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Abstract 


Energy resources and their utilization intimately relate to sustainable development. In attaining 
sustainable development, increasing the energy efficiencies of processes utilizing sustainable energy 
resources plays an important role. The utilization of renewable energy offers a wide range of 
exceptional benefits. There is also a link between exergy and sustainable development. A sustainable 
energy system may be regarded as a cost-efficient, reliable, and environmentally friendly energy 
system that effectively utilizes local resources and networks. Exergy analysis has been widely used in 
the design, simulation and performance evaluation of energy systems. 

The present study comprehensively reviews exergetic analysis and performance evaluation of a 
wide range of renewable energy resources (RERs) for the first time to the best of the author’s 
knowledge. In this regard, general relations (i.e., energy, exergy, entropy and exergy balance 
equations along with exergy efficiency, exergetic improvement potential rate and some thermo- 
dynamic parameters, such as fuel depletion ratio, relative irreversibility, productivity lack and 
exergetic factor) used in the analysis are presented first. Next, exergetically analyzed and evaluated 
RERs include (a) solar energy systems; (al) solar collector applications such as solar water heating 
systems, solar space heating and cooling, solar refrigeration, solar cookers, industrial process heat, 
solar desalination systems and solar thermal power plants), (a2) photovoltaics (PVs) and (a3) hybrid 
(PV/thermal) solar collectors, (b) wind energy systems, (c) geothermal energy systems, (cl) direct 
utilization (district heating, geothermal or ground-source heat pumps, greenhouses and drying) and 
(c2) indirect utilization (geothermal power plants), (d) biomass, (e) other renewable energy systems, 
and (f) country based RERs. Studies conducted on these RERs are then compared with the 
previously ones in tabulated forms, while the Grassmann (or exergy flow) diagrams, which are a very 
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useful representation of exergy flows and losses, for some RERs are given. Finally, the conclusions 
are presented. It is expected that this comprehensive study will be very beneficial to everyone involved 
or interested in the exergetic design, simulation, analysis and performance assessment of RERs. 

© 2006 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Achieving solution to environmental problems that we face today requires long-term 
potential actions for sustainable development. In this regard, renewable energy resources 
(RERs) appear to be the one of the most efficient and effective solutions [1]. 

RERs (i.e., solar, hydroelectric, biomass, wind, ocean and geothermal energy) are 
inexhaustible and offer many environmental benefits compared to conventional energy 
sources. Each type of renewable energy (RE) also has its own special advantages that make 
it uniquely suited to certain applications. Almost none of them release gaseous or liquid 
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Nomenclature 


City 
T 
< 


NS SSR ee by tO. ty SS. 
— 


Z 


amplitude, share of renewables in the total of each component 
area (m°) 

specific heat (kJ/kg K) 

coefficient of performance (dimensionless) 

energy rate (kW) 

specific exergy (kJ/kg) 

exergy rate (kW) 

exergetic factor, the sunlight dilution factor (dimensionless) 
exergy rate of fuel (kW) 

gravitational constant (m/s°) 

specific enthalpy (kJ/kg) 

heating value (kJ/kg) 

higher heating (gross calorific) value (kJ/kg) 

global irradiance (W/m?) 

rate of irreversibility, rate of exergy consumption (kW) 

rate of improvement potential (kW) 

enthalpy of phase change (kJ/kg) 

lower heating (net calorific) value (kJ/kg) 

mass flow rate (kg/s) 

number, index number (dimensionless) 

pressure (kPa) 

exergy rate of the product (kW) 

heat transfer rate (kW) 

renewable use by the residential-commercial sector in energy terms (kJ) 
ideal gas constant (kJ/kgK) 

specific entropy (kJ/kgK) 

entropy rate (kW/K) 

specific exergy index (dimensionless) 

period between local maxima and minima of the tidal record, time (s) 
temperature (°C or K) 

heat transfer coefficient (kW/m? K) 

speed (m/s) 

work (kJ) 

rate of work (or power) (kW) 

mol fraction (dimensionless) 

vertical distance from the water level of the reservoir to the reference height 
or average sea level (m) 

mass fraction (dimensionless) 


Greek letters 


Stefan—Boltzmann constant (W/m? K*) 
dilution factor (dimensionless) 
energy (first law) efficiency (dimensionless) 
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flow (specific) exergy (kJ/kg), maximum efficiency ratio (or exergy-to-energy 
ratio for radiation (dimensionless) 

interval 

density (kg/m?) 

proportionality constant (or quality factor or exergy coefficient) 

fuel depletion rate (dimensionless) 

exergy (second law) efficiency (dimensionless) 

productivity lack (dimensionless) 

relative irreversibility (dimensionless) 

specific humidity ratio (kgwater/k&air) 


air, actual 

absorber 

adsorbent 

ambient 

atmospheric 

available 

average 

beam 

Carnot 

cooking 
collector-evaporator 
chemical 

collector 

compressor 
condenser 

conversion 

cooker 

cooling 

natural direct discharge, diffuse 
drying air 

destroyed, destruction 
diffusive 

electrical, evaporator 
effective 

engine 

evaporator 

exergetic 

exergetic residential commercial 
fuel 

fan-coil 

vaporization 
generator 

geothermal district heating system 
generation, generated 
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ground heat exchanger 
heating, heat 

heat exchanger 

higher heating value 
heat pump 

heating value 
successive number of elements 
incompressible 

input 

inverter output 
location 

kinetic 

lower heating value 
maximum 

mechanical 

mixture 

overall 

overall electricity 
overall fuel 

overall residential 
output 

constant pressure, pump 
parabolic 

perfect 

physical 

potential 

power plant 
precipitation 

potential 

photovoltaic 

heat 

reinjected thermal water, refrigerant 
rational 

Rankine 
residential-commercial 
receiver 

reservoir 

residential overall 
solar 

solar collector 

space heating 

solar radiation 

system 

total, thermal 

thermal gradient 
transformation 

useful 
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v water vapor 

w water 

wa water/antifreeze solution 

wh wellhead, water heating 

wor without reinjected thermal water 
wr with reinjected thermal water 

0 dead (reference) state 

1 initial state 

2 final state 


pollutants during operation. In their technological development, the renewable ranges 
from technologies that are well established and mature to those that need further research 
and development [1,2]. 

Even though conventional sources, such as oil, natural gas and coal meet most 
of the energy demand at the moment, the role of RERs and their current advances have 
to take more relevance in order to contribute to energy supply and support the 
energy conservation (or efficiency) strategy by establishing energy management systems 
[3]. The use of RE offers a range of exceptional benefits, including: a decrease in 
external energy dependence; a boost to local and regional component manufacturing 
industries; promotion of regional engineering and consultancy services specializing in the 
utilization of RE; increased R&D, decrease in impact of electricity production and 
transformation; increase in the level of services for the rural population; creation of 
employment, etc. [4]. 

Dincer [5] reported the linkages between energy and exergy, exergy and the environment, 
energy and sustainable development, and energy policy making and exergy in detail. He 
provided the following key points to highlight the importance of the exergy and its 
essential utilization in numerous ways: (a) it is a primary tool in best addressing the impact 
of energy resource utilization on the environment. (b) It is an effective method using the 
conservation of mass and conservation of energy principles together with the second law of 
thermodynamics for the design and analysis of energy systems. (c) It is a suitable technique 
for furthering the goal of more efficient energy—resource use, for it enables the locations, 
types, and true magnitudes of wastes and losses to be determined. (d) It is an efficient 
technique revealing whether or not and by how much it is possible to design more efficient 
energy systems by reducing the inefficiencies in existing systems. (e) It is a key component 
in obtaining a sustainable development. 

Sustainable development does not make the world ‘ready’ for the future generations, but 
it establishes a basis on which the future world can be built. A sustainable energy system 
may be regarded as a cost-efficient, reliable, and environmentally friendly energy system 
that effectively utilizes local resources and networks. It is not ‘slow and inert’ like a 
conventional energy system, but it is flexible in terms of new techno-economic and political 
solutions. The introduction of new solutions is also actively promoted [6]. 

An exergy analysis (or second law analysis) has proven to be a powerful tool in the 
simulation thermodynamic analyses of energy systems. In other words, it has been widely 
used in the design, simulation and performance evaluation of energy systems. Exergy 
analysis method is employed to detect and to evaluate quantitatively the causes of the 
thermodynamic imperfection of the process under consideration. It can, therefore, indicate 
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the possibilities of thermodynamic improvement of the process under consideration, but 
only an economic analysis can decide the expediency of a possible improvement [7,8]. 

The concepts of exergy, available energy, and availability are essentially similar. The 
concepts of exergy destruction, exergy consumption, irreversibility, and lost work are also 
essentially similar. Exergy is also a measure of the maximum useful work that can be done 
by a system interacting with an environment which is at a constant pressure Py and a 
temperature To. The simplest case to consider is that of a reservoir with heat source of 
infinite capacity and invariable temperature Tọ. It has been considered that maximum 
efficiency of heat withdrawal from a reservoir that can be converted into work is the 
Carnot efficiency [9,10]. 

Although numerous studies have been conducted on the energetic analysis and 
performance evaluation of RERs by using energy analysis method in the literature, very 
limited review papers have appeared on exergy analysis and performance assessment of 
RERs. In this regard, Koroneos et al. [11] dealt with the three kinds of RERs in terms of 
exergetic aspects, namely (i) exergetic analysis of a solar thermal power system by 
presenting a study conducted by Singh et al. [12], (ii) exergy analysis of geothermal power 
systems by presenting and evaluating the Larderello—Farinello—Valle Secolo Geothermal 
Area (Tuscany, Italy) studied by Bettagli and Bidini [13], and (iii) exergy analysis of wind 
energy systems by giving the relations used in the analysis and evaluating a wind energy 
system. They also made a comparison between renewable and non-RE sources, and 
concluded that some of the systems appear to had high efficiencies, and in some cases they 
are greater than the efficiency of systems using non-RE sources. In other cases, like the 
conversion of solar energy to electricity, the efficiencies were lower, in order to meet the 
electricity needs of cities. Hermann [14] identified the primary exergy reservoirs that supply 
exergy to the biosphere and quantified the intensive and extensive exergy of their derivative 
secondary reservoirs, or resources. Exergy relations of cosmic radiation exchange, wind, 
ocean surface waves, precipitation, ocean thermal gradient, tides and geothermal were also 
briefly presented. 

The present study differs from the previous ones due to the facts that: (1) This covers a 
comprehensive exergetic analysis and performance evaluation of RERs. (ii) This includes a 
wide range of RERs such as solar, wind, geothermal and biomass along with their 
subsections (i.e., for geothermal, its direct and indirect applications) as well as hybrid 
systems. (iii) This presents a comparison between previously conducted studies in tabulated 
forms, and (iv) This summarizes exergetic utilization efficiency of RERs as done by 
Hepbasli and Utlu [15] for the first time for Turkey. In this regard, the structure of the 
paper is organized as follows: The first section includes the introduction; Section 2 deals 
with energy and exergy modeling, while general relations are described in Section 3; 
exergetic analysis and evaluating the RERs are treated in Section 4 in more detail by 
applying the general relations to various RERs; and the last section concludes. 


2. Energy and exergy modeling 


Dincer et al. [16] reported that, to provide an efficient and effective use of fuels, it is 
essential to consider the quality and quantity of the energy used to achieve a given 
objective. In this regard, the first law of thermodynamics deals with the quantity of energy 
and asserts that energy cannot be created or destroyed, whereas the second law of 
thermodynamics deals with the quality of energy, i.e., it is concerned with the quality of 
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energy to cause change, degradation of energy during a process, entropy generation and 
the lost opportunities to do work. More specifically, the first law of thermodynamics is 
concerned only with the magnitude of energy with no regard to its quality; on the other 
hand, the second law of thermodynamics asserts that energy has quality as well as 
quantity. By quality, it means the ability or work potential of a certain energy source 
having certain amount of energy to cause change, i.e., the amount of energy which can be 
extracted as useful work which is termed as exergy. First and second law efficiencies are 
often called energy and exergy efficiencies, respectively. It is expected that exergy 
efficiencies are usually lower than the energy efficiencies, because the irreversibilities of the 
process destroy some of the input exergy. 

Exergy is the expression for loss of available energy due to the creation of entropy in 
irreversible systems or processes. The exergy loss in a system or component is determined 
by multiplying the absolute temperature of the surroundings by the entropy increase. 
Entropy is the ratio of the heat absorbed by a substance to the absolute temperature at 
which it was added. While energy is conserved, exergy is accumulated [17]. 

Exergy analysis provides a method to evaluate the maximum work extractable from a 
substance relative to a reference state (i.e., dead state). This reference state is arbitrary, but 
for terrestrial energy conversion the concept of exergy is most effective if it is chosen to 
reflect the environment on the surface of the Earth. The various forms of exergy are due to 
random thermal motion, kinetic energy, potential energy associated with a restoring force, 
or the concentration of species relative to a reference state. In order to establish how much 
work potential a resource contains, it is necessary to compare it against a state defined to 
have zero work potential. An equilibrium environment which cannot undergo an energy 
conversion process to produce work is the technically correct candidate for a reference 
state [14]. 

It should be noticed that exergy is always evaluated with respect to a reference 
environment (i.e. dead state). When a system is in equilibrium with the environment, the 
state of the system is called the dead state due to the fact that the exergy is zero. At the 
dead state, the conditions of mechanical, thermal, and chemical equilibrium between the 
system and the environment are satisfied: the pressure, temperature, and chemical 
potentials of the system equal those of the environment, respectively. In addition, the 
system has no motion or elevation relative to coordinates in the environment. Under these 
conditions, there is neither possibility of a spontaneous change within the system or the 
environment nor an interaction between them. The value of exergy is zero. Another type of 
equilibrium between the system and environment can be identified. This is a restricted form 
of equilibrium, where only the conditions of mechanical and thermal equilibrium (thermo- 
mechanical equilibrium) must be satisfied. Such state is called the restricted dead state. At 
the restricted dead state, the fixed quantity of matter under consideration is imagined to be 
sealed in an envelope impervious to mass flow, at zero velocity and elevation relative to 
coordinates in the environment, and at the temperature To and pressure Po taken often as 
25°C and | atm [18]. 


3. General relations 
For a general steady state, steady-flow process, the four balance equations (mass, 


energy, entropy and exergy) are applied to find the work and heat interactions, the rate of 
exergy decrease, the rate of irreversibility, the energy and exergy efficiencies [16,19—21]. 
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3.1. Mass, energy, entropy and exergy balances 


The mass balance equation can be expressed in the rate form as 


So iin = So touts (1) 


where 7i2 is the mass flow rate, and the subscript in stands for inlet and out for outlet. 
The general energy balance can be expressed below as the total energy inputs equal to 
total energy outputs. 


X Èn = XO Eou- (2) 


In the absence of electricity, magnetism, surface tension and nuclear reaction, the 
total exergy of a system Ex can be divided into four components, namely (i) physical 


exergy Èx", (ii) kinetic exergy Èx“^, (iii) potential exergy Èx”, and (iv) chemical exergy 
- CH 
Ex [16]. 
Ex = Èx" + Ex" + Èx" + Èx". (3) 
Although exergy is extensive property, it is often convenient to work with it on a unit of 
mass or molar basis. The total specific exergy on a mass basis may be written as follows: 
ex = exPH + exKN + ex’! 4 ex™, (4) 


The general exergy balance can be written as follows: 


> EXin = 5 Èxout = D Exdest (Sa) 


or 
EXheat a Exwork F EX mass,in = EX mass,out = Ex¢est (Sb) 
with 
A To\ ; 
EXneat =% (1 = D Qr» (6a) 
k 
Exwork = W, (6b) 
EXimass,in = Erioni out> (6c) 
Exmass,out = ZMoutW out> (6d) 


where Q, is the heat transfer rate through the boundary at temperature T+ at location k 


and W is the work rate. 
The flow (specific) exergy is calculated as follows: 


y = (h — ho) — Tos — so), (7) 


where h is enthalpy, s is entropy, and the subscript zero indicates properties at the 
restricted dead state of Pp and Tp. 
The rate form of the entropy balance can be expressed as 


Sin = Sout + Sx = 0, (8) 
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where the rates of entropy transfer by heat transferred at a rate of Q, and mass flowing at a 
rate of m are Sheat = Ò; /Tk and S mass = MS, respectively. 

Taking the positive direction of heat transfer to be to the system, the rate form of the 
general entropy relation given in Eq. (8) can be rearranged to give 


: Ò; 
Sgen = ZMoutSout = LMinSin a hn: (9) 

k 
Also, it is usually more convenient to find Sien first and then to evaluate the exergy 


destroyed or the irreversibility rate J directly from the following equation, which is called 
Gouy-Stodola relation [22]: 


fate Tosa (10) 


The specific exergy (flow exergy) of an incompressible substance (i.e., water) is given 
by [22] 


0 


by =C(T-To-Toln 7). (i) 
The total flow exergy of air is calculated from [23] 


Wat = (Cpa + @Cpy)Tol(T/To) — 1 — In(T/To)] + (1 + 1.60780)Ra To In(P/Po) 


(1 + 1.60780) In[(1 + 1.607899) 
R,T : 12 
= ofa + 1.6078@)] + 1.60780 In(w/coo) 2) 
where the specific humidity ratio is 
œ = my, /Ma. (13) 


Assuming air to be a perfect gas, the specific physical exergy of air is calculated by the 
following relation [24] 


T P 
Waper = Cra( T- To — To In 7) + RaTo In Pa (14) 


3.2. Energy and exergy efficiencies 


Numerous ways of formulating exergetic (or exergy or second-law) efficiency 
(effectiveness, or rational efficiency) for various energy systems are given in detail 
elsewhere [20]. It is very useful to define efficiencies based on exergy (sometimes called 
Second Law efficiencies). Whereas there is no standard set of definitions in the literature, 
two different approaches are generally used—one is called 1 “brute-force”, while the other 
is called “functional.” [25]. 


e A “brute-force” exergy efficiency for any system is defined as the ratio of the sum of all 
output exergy terms to the sum of all input exergy terms. 

e A “functional” exergy efficiency for any system is defined as the ratio of the exergy 
associated with the desired energy output to the exergy associated with the energy 
expended to achieve the desired output. 
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Here, in a similar way, exergy efficiency is defined as the ratio of total exergy output to 
total exergy input, i.e. 


= EXoutput =1 EXdest 


(15) 


= > 
EXinp ut Exinput 


where “output or out” stands for “net output” or “product” or “desired value” or 
“benefit”, and “input or in” stands for “given” or “used” or “fuel”. 

It is clear that the brute-force definition can be applied in a straightforward manner, 
irrespective of the nature of the component, once all exergy flows have been determined. 
The functional definition, however, requires judgment and a clear understanding of the 
purpose of the system under consideration before the working equation for the efficiency 
can be formulated [25]. 


3.3. Exergetic improvement potential 


Van Gool [26] has also proposed that maximum improvement in the exergy 
efficiency for a process or system is obviously achieved when the exergy loss or 
irreversibility (Ein — Eou) is minimized. Consequently, he suggested that it is useful 
to employ the concept of an exergetic “improvement potential” when analyzing different 
processes or sectors of the economy. This improvement potential in the rate form, denoted 
IP, is given by 


IP = (1 — 8)(Exin — EXout). (16) 


3.4. Some thermodynamic parameters 


Thermodynamics analysis of RERs may also be performed using the following 
parameters [27]: 
Fuel depletion ratio: 


rf 
j=. (17) 
Fr 
Relative irreversibility: 
I 
=>. 18 
=i (18) 
Productivity lack: 
fi 
==. 19 
č P, (19) 
Exergetic factor: 
F; 
Íi (20) 
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4. Exergetic analysis and evaluation of renewable energy resources 
4.1. Solar energy systems 


In evaluating the performance of solar energy systems using exergy analysis method, 
calculation of the exergy of radiation is very crucial. However, its calculation is a problem 
of unquestionable interest, since exergy represents the maximum quantity of work that can 
be produced in some given environment (usually the terrestrial environment, considered as 
an infinite heat source or sink) [28]. Over a period of more than 20 years, many papers 
including various approaches to this calculation have been published, as reported in 
[29,30]. Among these, the first one [31], which is called by Petela [29,32] the maximum 
efficiency ratio (or exergy-to-energy ratio for radiation) for determining an exergy of 
thermal emission at temperature T, is as follows: 


1/T\f ATs 
Vadis =1+ 3 (2) > Fe (21) 


where T was taken to equal the solar radiation temperature (7,) with 6000 K in exergetic 
evaluation of a solar cylindrical-parabolic cooker and a solar parabolic-cooker by Petela 
[32] and Ozturk [33], respectively. Eq. (21) was also derived by Szargut [34] who presented 
a simple scheme of a reversible prime mover (different from that of Petela) transforming 
the energy of radiation into mechanical or electrical work and took into account that the 
solar radiation had a composition similar to that of a black body. 

The history of Eq. (21) is as follows [35]: It was first derived by Petela [31] and 
independently by Landsberg and Mallinson [36] and Press [37], while Spanner [38] 
proposed Petela’s equation for the direct solar radiation exergy [29]. As a consequence, it is 
sometimes called the Landsberg efficiency [34], the Petela-Landsberg upper bound [39], or, 
more appropriately, the Petela—Press—Landsberg and Mallinson formula [28]. 

Petela [29] reported that in the existing literature in the above field, the most discussed 
formulae for heat radiation exergy were those derived by him [31], Spanner [38] and Jeter 
[40] given in Eqs. (21), (24b), and (24c), respectively. As for as the application of these 
equations to various thermal systems by some investigators, Chaturvedi et al. [41] used 
Eq. (21) in the performance assessment of a variable capacity direct expansion solar- 
assisted heat pump system for domestic hot water application. Eskin [42] considered 
Eq. (24b) in evaluating the performance of a solar process heating system, while Singh et 
al. [12], and Ucar and Inalli [43] considered Eq. (24c) in the exergetic analysis of a solar 
thermal power plant and the exergoeconomic analysis of a solar-assisted heating system, 
respectively. 

Petela [29] also listed three various limiting energy efficiencies of utilization of the 
radiation matter in Table 1 [29] and explained them as follows: “Bejan’s conclusion was 
that all theories (Spanner’s, Jeter’s and Petela’s) concerning the ideal conversion of thermal 
radiation into work, although obtaining different results, were correct. Although really 
correct, they had incomparable significance for the true estimation of the radiation exergy 
value. All of the efficiencies assumed work as an output. However, Petela’s work was equal 
to the radiation exergy, Jeter’s work was the heat engine cycle work and Spanner’s work 
was an absolute work. As an input, Petela and Spanner assumed the radiation energy, 
whereas Jeter assumed heat [29].” 
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Table 1 
Numerators (output) and denominators (input) of the limiting energy efficiency of radiation utilization by three 
different researchers [29] 


Investigators Input Output Unified efficiency expression 
Petela [31] Radiation energy Useful work = radiation exergy 1/T)\* 4To 
l+5(—) -z7 QI’) 

3\T 3T 

Spanner [38] Radiation energy Absolute work 1- 4To (24b) 
3T 

Jeter [40] Heat Net work of a heat engine i= To (24c) 
T 


Svirezhev et al. [44] applied a new concept regarding the exergy of solar radiation to the 
analysis of satellite data describing the seasonal (monthly) dynamics of four components of 
the global radiation balance: incoming and outgoing short and long-wave solar radiation. 
They constructed maps of annual mean exergy and some derivative values using real 
satellite data. They concluded that exergy could be a good indicator for crucial regions of 
the oceans, which were indicated where its values were maximal. 


4.1.1. Solar collector applications 
The instantaneous exergy efficiency of the solar collector can be defined as the ratio of 
the increased water exergy to the exergy of the solar radiation [33]. In other words, it is a 
ratio of the useful exergy delivered to the exergy absorbed by the solar collector [12]. 
Exu 


é =-= 22 
scol EX.col ( ) 
with 
Ex, = My [hy out = hwin) = To(Sw,out = Sw,in)] (23a) 
or 
T 
Exu = my Cy Poo = Twin) — To (in rat) > (23b) 
w,in 
: : T Tw ar) 
Ex, = 1 Ln— 23c 
5 Qu | T wont ag Twin ( Twin ( ) 
and 
ExXgcol = Al pW srad.maxo (24a) 


where the subscripts “in” and “out”? denote the inlet and outlet of the solar collector, 
respectively, while W.-admax May also be calculated using Eqs. (21), (24b) and (24c) given in 
Table 1 [29]. 

The exergy of solar radiation with beam (Ipe) and diffuse (74) components for parabolic 
collectors is also given by the following equation [45], which was used by Eskin [46] in the 
performance evaluation of a cylindrical parabolic concentrating collector. 


4T 4T 
W srad,be,d = Ibe (1 = T) +a (: = T) (25) 
s 


606 A. Hepbasli / Renewable and Sustainable Energy Reviews 12 (2008) 593-661 


Ts 
ro 0.9562 + 0.2777 In(1/K«) + 0.0511, (26) 


ok 
S 


where x is the dilution factor of diffuse radiation [47]. 


4.1.1.1. Solar water heating systems. Solar water heater is the most popular means of 
solar energy utilization because of technological feasibility and economic attraction 
compared with other kinds of solar energy utilization. Fig. 1 illustrates a schematic of a 
typical domestic-scale solar water heater, which mainly consists of a solar collector and a 
storage barrel [48]. 

Exergy efficiency of the domestic-scale solar water heater may be calculated using 
Eq. (15) as follows: 


ÉXoutput (27) 


Es heater = = > 
Í EXsun 


where exergy from the storage barrel to the end-user (output exergy) may be written as 
follows [48], assuming the temperature distribution in the storage barrel is linear: 


Top +T T 
Bsouput = ity Co ( top 2 bottom To) iy CoTo( m21) 


T bottom Tomw Cw In Top 
Top = T'bottom T bottom 


(28) 


and Exgun is the exergy from sun (input exergy) and may be found using Eq. (24a). 
The collector exergy efficiency is calculated by the following equation, while it may be 
found using Eqs. (22), (23b) and (24). 


EXcol,tank 
EXsun 


> (29) 


Esco] = 


where EXcoltank is exergy from the collector to the storage barrel and may be written in a 
similar way given by Eq. (23b) as follows [48]. 


Twou 
EXcol,tank = My Cy (Too = To) — To (1m Troe) : (30) 
0 


ave 


Storage 
barrel 


— 


Fig. 1. A typical domestic-scale solar water heater [48]. 
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Xiaowu and Ben [48] performed an exergy analysis of a domestic-scale solar heater. They 
used the data obtained from Domestic-Scale Solar Water Heaters Performance Testing 
Station, located in Yun Nan province, China as follows: Top and bottom hot water 
temperatures of 50 and 38 °C in the storage barrel, respectively, 2.5m? of collector area, an 
average ambient temperature of 25°C, an average ambient water temperature of 24°C, 
466 W/m? of whole-day average solar radiation and 196.4kg of water storing log capacity 
for the storage barrel. They calculated the exergy efficiency value of the solar water heater 
using the above data and found to be 0.77%. The main conclusions drawn from the results 
of their study were as follows: The exergy efficiency of domestic-scale water heater was small 
as the output energy was of low quality. The storage barrel had large exergy losses. It was 
worthwhile to design a new style of storage barrel for reducing unnecessary mixing of water. 
In order to improve the exergy efficiency of domestic-scale water heater, a judicious choice 
of width of plate and layer number of cover was necessary [48]. 


4.1.1.2. Solar space heating and cooling 

4.1.1.2.1. Air systems. Air systems are indirect water heating systems that circulate air 
via ductwork through the collectors to an air-to-liquid heat exchanger. In the heat 
exchanger, heat is transferred to the potable water, which is also circulated through the 
heat exchanger and returned to the storage tank. A double storage tank system is shown in 
Fig. 2 [48]. This type of system is used most often, because air systems are generally used 
for preheating domestic hot water and thus auxiliary is used only in one tank. The main 
advantage of the system is that air does not need to be protected from freezing or boiling, is 
noncorrosive, and is free. The disadvantages are that air handling equipment (ducts and 
fans) need more space than piping and pumps, air leaks are difficult to detect, and parasitic 
power consumption is generally higher than that of liquid systems [49]. 


AIR-TO-WATER 
OUT DOOR HEAT EXCHANGER 


AUXILIARY 
STORAGE (PREHEAT) 
TANK STORAGE 
TANK 


COLD WATER 


Fig. 2. A schematic of air system with a double storage tank drawn by Xiaowu and Ben [48]. 
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The relations used to evaluate air systems using exergy analysis method are presented by 
the author since any exergetic assessment based on the actual data could not be found. The 
main parts to be included in the exergy analysis are (J) air solar collector, (II) air-to-water 
hat exchanger, (III) fan, (IV) storage tank and (V) circulating pump. The exergy 
destruction and exergy efficiency relations for each component of the whole air system are 
given below using Eqs. (5a) and (15), respectively. 


I. Air solar collector: 


Exdest,1 = EXin air = Exoutair + Exs, (31) 


EXout air — EXin air 
a= = ; (32) 
EXscol 


where EXin air and ÈXoutair are calculated using Eqs. (6c), (6d), (12) and (14), while 
Exgco is found from Eq. (24a). 
II. Air-to-heat exchanger: 


EXaest 11 = (Ežnot,air,in ag EXcold,w,in) a (EXhotair,out + EXcold,w,out)» (33) 


ExXcold,w,out ca EXcold,w.in = Meoldw(Weold,w,out Weold,w.in) (34) 


=z g =e : 
EXhotair,in = EXhotair,out Mpotair Wpotair,in = Whotainout) 


III. Fan: 
EX¢est,111 = Wane F Mair(Woutair = Win air)> (35) 
Exout — Exin  Mair(Woutair — Vinai 
én = Xout Xin = air(Wo tair Win air) (36) 
W tan W tan 
with 
; y2. 
W fan = Mair (Aoutair = hin air) + n . (37) 
IV. Storage tank: 
Exzest,v and ery may be calculated in a similar way given by Eqs. (33) and (34). 
V. Circulating pump: 
ExXdest,vi = Mw(Win — Wouta) + W purge (38) 
ec EXouta om Exin = Mw (outa = Win) (39) 


W pump. W pump. 
with 


W pump.e = Wcities Miata ts) (40) 
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where Npump,e and Ypump,mech are the circulating pump motor electrical and circulating 
pump mechanical efficiencies. 


The overall exergy destruction in the whole air system is calculated as follows, while the 
overall exergy efficiency is found on a product/fuel basis given by Eq. (15). 


Exdest,o = ExXdest,I + EXdest,11 + ExXdest,IIl + ExXdest,IV + ExXdest,v- (41) 


4.1.1.2.2. Water systems. There are many variations of systems used for solar space 
heating. The exergetic analysis of these systems may be performed using relations given in 
Sections 4.1.1.1 and 4.1.1.2.1. 

4.1.1.2.3. Heat pumps. Heat pumps are preferred and widely used in many 
applications due to their high utilization efficiencies compared to conventional heating 
and cooling systems. There are two common types of heat pumps: air-source heat pumps 
and ground-source (geothermal) heat pumps (GSHPs). The main components of these 
systems are a compressor, a condenser, an expansion valve and an evaporator. In addition, 
circulating pumps and fans are used. GSHPs are exergetically studied in more detail in 
Section 4.3.2.2 and relevant exergy relations are presented in a tabulated form, while 
studies conducted by various investigators are compared. The relations given in these 
sections may be used in modeling and exergetic evaluation for heat pump systems. In case 
that a solar assisted heat pump with direct expansion of the refrigerant within the solar 
collector is used, the collector-evaporator part may be analyzed using Eq. (24a) for solar 
radiation input. 

Tsaros et al. [50] studied on exergy analysis of heat pumps. Exergy transports between 
the components, to the loads and to the environment were assessed. The exergy 
consumption in each component were also determined. Typical results were presented for 
the operation of a nominal 3.5 ton residential air-to-air heat pump, in both the heating 
and cooling modes. Exergy efficiency values for various system components are listed in 
Table 2 [50]. Of the major components, the outdoor coil had the highest exergy efficiency, 
followed by the compressor and the indoor coil. The exergy efficiency of the indoor coil 
varied greatly with load and exceeded that of the compressor under some conditions. The 
fans, with the values less than 5%, were the least efficient [50]. 

Bilgen and Takahashi [51] carried out exergy analysis of heat pump—air conditioner 
systems. The experimental system, shown in Fig. 3 [51], was Matsushita room air 
conditioner CS-XG28M. The irreversibilities due to heat transfer and friction were 
considered. The coefficient of performance based on the first law of thermodynamics as a 
function of various parameters, their optimum values, and the efficiency and coefficient of 
performance based on exergy analysis were derived. Based on the exergy analysis, a 
simulation program was developed to simulate and evaluate experimental systems. It was 
found that the exergy efficiency varied from 37% to 25% both a decreasing function of 
heating or cooling load. The exergy destructions in various components were also 
determined for further study and improvement of its performance. It was suggested that to 
improve the performance of the heat pump system, each component could be further 
studied from exergy usage and economics points of view. Eventually, a system exergy- 
economics optimization study could be carried out to identify each component’s relative 
importance with respect to operational conditions. 
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Table 2 
Exergy efficiency values of the system components [50] 


Component Exergy efficiencies in % at various outdoor temperatures 
—18.75°C (17°F) 8.33 °C (47°F) 35°C (95°F) 
Indoor coil 82.6 57.0 35.3 
Compressor, 4-way valve, outdoor coil, throttle 36.0 40.8 29.8 
Compressor, 4-way valve, outdoor coil 44.0 47.7 35.2 
Throttle 81.6 85.5 84.5 
Compressor, 4-way valve 47.9 55.0 SRA 
Outdoor coil 92.0 86.9 65.6 
Compressor 50.5 58.5 58.3 
4-way valve 94.7 93.9 92.1 
Outdoor fan and motor 4.41 4.41 4.41 
Indoor fan and motor 3.49 3.49 5.81 
Overall without fans 30.1 23.4 10.7 
Overall 21.7 18.1 8.43 


Compressor includes associated hardware, while 4-way valve includes piping to coils. 
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Fig. 3. Schematic of the experimental set up showing measurement system and instrumentation drawn by Bilgen 
and Takahashi [51]. 


Aprea et al. [52] studied on the performances of a vapor compression experimental 
system working both as a water chiller and as a heat pump, using as refrigerant fluids R22 
and its substitute R417A (R125/R134a/R600 46.6/50/3.4% in mass). This kind of system 
has been used in industrial processing or air conditioning plants where a supply of 
refrigerated and heated water was required. The experimental system was principally made 
up of a unit that could operate both as refrigeration system and heat pump. It was made up 
of a hermetic compressor Scroll, a plate-type water heat exchanger inserted in a water tank, 
a finned tube air heat exchanger, two thermostatic expansion valves. 
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Aprea et al. [52] evaluated the efficiency defect for each component of the system, 
considering the ratio between the exergy destroyed in each one and the exergy required to 
sustain the process (i.e., the electrical power supplied to the compressor). The efficiency 
defects of R417A were always higher than the defects of R22 of about 6% related to the 
valve and the heat exchangers, and of about 10% referring to the compressor. Moreover, it 
was experimentally verified that the COP of R22 was higher than that of R417A also of 
about 18% related to the water chiller system and of about 15% referring to the heat 
pump. 

Ma and Li [53] exergetically analyzed a heat pump system with economizer coupled with 
scroll compressor illustrated in Fig. 4 [53] and derived expressions for exergy losses of the 
processes in the system, exergy loss ratio and efficiency of the system. They also determined 
the exergetical losses and efficiency of this heat pump system using experimental data, 
while they assessed characteristics of exergy variation and transport and exergetic 
efficiency of the system. The experimental system studied is explained in more detail 
elsewhere [54]. 

Using Eq. (15), Ma and Li [53] expressed exergy efficiency as a ratio of the effective 
exergy output of the heat pump (equal to the heat load exergy in the condenser plus the 
cool capacity exergy in the evaporator) to the exergy of the electrical energy as follows: 

T EXeff, output = Exp + EXcool 


e= 42 
EXinput Ex, i ) 


From the distribution of the exergy loss, it was found that the maximum loss occurred in 
the compressor, accounting for about 82% of the total exergy destruction. This value 
would increase when the temperature difference between the condensation and the 
evaporation could be increased. The compressor was a bottleneck for the efficiency of the 
heat pump system, and increasing energy efficiency ratio for the compressor was crucial to 
the improvement of the heat pump system. For hermetic compressor, there were two 
mechanisms for exergy to consumption: one is was the process in which electrical power 
was converted to mechanical work in motor, and the other was the process of the 
mechanical work being converted to the pressure energy of the refrigerant in the 
compressor. The exergy losses in the evaporator and condenser were found to be 8.9% and 


scroll compressor 
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expansion valve A 


condenser 


Fig. 4. Schematic of a heat pump system with economizer coupled with scroll compressor drawn by Ma 
and Li [53]. 
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5.3% of the total loss, respectively, and were 10.8% and 6.5% of the loss in the 
compressor. Although they were relatively small compared with the loss in the compressor, 
minimizing the exergy losses in evaporator and condenser was also important for the 
improvement of the efficiency of the heat pump, and improving the evaporator was more 
beneficial for the overall system efficiency than improving the condenser. The exergy losses 
existing in the essential components of the heat pump system such as the compressor, the 
evaporator and the condenser were predominant in the total exergy loss, namely 96.2% of 
the total. In summary, the best gain in the heat pump system performance could be 
obtained through improvement in the compressor, followed by the evaporator and 
condenser. Using Eq. (42), exergy efficiency of the air-source heat pump system was 
determined to vary approximately from 33% to 42% at the evaporation temperatures of 
—25 to —12°C with a fixed condensing temperature of 45°C [53]. 

Badescu [55] proposed a model for the heating system of an ecological building 
whose main energy source was solar radiation. The most important component of the 
heating system was a vapor compression heat pump. The building’s electric energy was 
provided by a photovoltaic (PV) array or from the AC grid, when necessary. Exergy 
analysis of the heat pump was performed. The state parameters and the process quantities 
were evaluated by using, as input, the building thermal load. It was reported that most of 
the exergy losses occurred during compression and condensation. Preliminary results 
indicated that the PV array could provide all the energy required to drive the heat pump 
compressor, if an appropriate electrical energy storage system was provided. In the 
analysis, meteorological data measured by the Romanian Meteorological and Hydro- 
logical Institute during 1961 in Bucharest (latitude 44.5°N, longitude 26.2°E, altitude 
131 m above sea level) were used, while the exergy efficiency of the heat pump, was found 
using the following equation: 


fee ee (43) 
Weomp,a 
where Weomp,a and Weomp,c are the specific mechanical works entering the actual and 
reversible (Carnot) compressor, respectively. Using Eq. (43), exergy efficiency of the heat 
pump was found to be 15.53% and 18.37% for the refrigerants of R114 and R12, 
respectively, based on the results obtained in the cold day of January 22, 1961 at 6.00 a.m. 
in Bucharest, when the ambient temperature was 255.7 K [55]. 

Cervantes and Torres-Reyes [56] performed an experimental study of a solar assisted 
heat pump with direct expansion of the refrigerant within the solar collector and evaluated 
the performance of this system using exergy analysis method. The condensation of the 
refrigerant took place in an R22/air condenser made out of stainless-steel, with a heat 
transfer area of 9.8 m”. The compression stage was carried out by a hermetic compressor, 
and the expansion of R-22 occurred in a thermal expansion valve. The refrigerant was 
evaporated inside an uncovered solar flat collector, with a collection area of 364.5 m”. 
Experimental tests were performed at a solar radiation range of 200-1100 W/m? and 
ambient temperatures varying from 20 to 32°C. The actual amount of work supplied to 
the compressor, given by the electricity consumption, ranged from 1.1 to 1.36kW. The 
maximum exergy efficiency, defined as the ratio of the outlet to the inlet exergy flow 
in every component of the heat pump cycle, was determined taking into account the 
typical parameters and performance coefficients. The results of this exergy analysis 
pointed out that the main source of irreversibility could be found in the evaporator of 
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the heat pump (ie.,the solar collector) emphasizing that incoming solar radiation, as 
shown in Fig. 5 [56]. A full analysis of the inlet and outlet exergy flows and the main 
parameters associated to them, for each component of the solar assisted heat pump, is 
given in Ref. [57]. 

Ucar and Inalli [43] developed an exergoeconomic model for analysis and optimization 
of solar-heating systems with residential buildings. The model system consisted of flat plate 
solar collectors, a heat pump, an underground storage tank and a heating load, as shown 
in Fig. 6 [43]. The heat produced by the solar collectors was stored in underground storage 
tanks throughout the year. This heat was only extracted during the heating season. It was 
found from the exergetic analysis results for the cylindrical tank system that the largest 
exergy losses occurred in the house at 78.7%, followed by the heat pump at 40.9%, the 
collector at 39.7% and the storage at 19.8%. 


Collector - evaporator Tce 


W, Compressor 


Fig. 5. Grassmann (exergy flow) diagram of the solar assisted heat pump system with direct expansion of the 
refrigerant within the solar collector drawn by Cervantes et al. [56]. 
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Fig. 6. Schematic diagram of a solar-assisted heating system for residential buildings drawn by Ucar and Inalli 
[43]. 
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4.1.1.3. Solar refrigeration 

4.1.1.3.1. Solar-driven ejector refrigeration systems. A solar-driven ejector refrigera- 
tion system consists of a solar collector subsystem and an ejector refrigeration subsystem. 
The major components in the refrigeration cycle are an ejector, a condenser, a generator, 
an evaporator, an expansion device and a pump. Fig. 7 [58] illustrates a schematic of this 
system, which was energetically analyzed by Pridasawas and Lundqvist [58] in more detail. 
In assessing this system using exergy analysis method, the relations given in Section 4.3.2.2 
and additionally the exergetic efficiency definitions explained below may be used. 

The exergy inputs and outputs may take different forms, e.g., radiation, heat or 
electricity. Thus, the exergetic efficiency of the ejector refrigeration cycle is defined as the 
ratio of the cooling exergy (at the evaporator) and the exergy inputs to the generator and 
the pump: 


EXey; 
Ej = a (44) 
Ex, + Expe 
with 
a iy 
Ex, = @,(1 = D] (45) 
g 


where Ò, is energy input to the generator and T, is the generating temperature [58]. 
The system exergetic efficiency is defined as the ratio of the cooling exergy and the 
exergy of the available solar heat input plus the electricity input to the pump: 


EXevap 
=e 46 
Esys EX<h + Expe ( ) 
with 
. . To 
Exsh = Q ava (1 — T ), (47) 
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Fig. 7. Ejector refrigeration cycle drawn by Pridasawas and Lundqvist [58]. 
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where Q,,,, is available heat for the process in the solar collector and Tyco is the average 
temperature of the solar collector [58]. 

Pridasawas and Lundqvist [58] used the following relation established by De Vos [59] for 
calculating the exergy input to the collector per unit area as a function of the exergy 
emitted from the sun minus the albedo of the earth and the radiation emitted from the 
solar collector: 

EXsoot = foT? +(1—f)6T hance — 9T s 


scol? 


(48) 


where f is the sunlight dilution factor, which equal to 2.16 x 107° on the earth, ø is 
Stefan—Boltzmann constant (5.67 x 1078 W/m’ K*), T, is the temperature of the sun (K) 
and Tyianet is the temperature of the planet (K). 

The exergy loss during the transformation from solar radiation to heat on the solar 
collector is given as [58]: 


I scol,trans = Exs = Exs hn. (49) 


The exergy loss from the exergy input to the solar collector to the working fluid can be 
calculated using the following equations [58]: 


I scol = Exsh — Exy (50) 
with 
Ex, = @,(1 aon | (51) 
T scol 


Pridasawas and Lundqvist [58] evaluated the performance of an ejector refrigeration 
cycle driven by solar energy using exergy analysis method. The analysis was based on the 
following conditions: a solar radiation of 700 W/m’, an evaporator temperature of 10°C, a 
cooling capacity of 5kW, butane as the refrigerant in the refrigeration cycle and ambient 
temperature of 30°C as the reference temperature. Irreversibilities occurred among 
components and depend on the operating temperatures. The most significant losses in the 
system were in the solar collector and the ejector. The latter decreased inversely 
proportional to the evaporation temperature and dominated the total losses within the 
system. For the above operating conditions, the optimum generating temperature was 
about 80 °C. The exergy losses in the ejector refrigeration cycle is shown in Fig. 8 [58]. The 
highest loss with about 50.5% was found in the solar collector, followed by the losses in the 
ejector, which was 15.8% and in the condenser with 10.6% (1.49kW = 0.89 kW loss and 
0.60kW rejected). On the component basis of the ejector refrigeration subsystem, the 
largest loss (32%) occurred in the ejector followed by the generator (22%), the condenser 
(21%). The rest of the exergy loss occurred in the pump, the evaporator and the expansion 
valve. The exergetic efficiency values for the solar collector and refrigeration cycle were 
14.53 and 4.39%, respectively, while the overall exergetic efficiency was 0.66%. 

4.1.1.3.2. Absorption units. Sozen et al. [60] investigated a prototype of an 
aqua—ammonia absorption heat pump system using solar energy. Fig. 9 [60] illustrates a 
schmatic of this system. The aqua-ammonia absorption heat pump consisted of condenser, 
evaporator, absorber, generator, pump, expansion valves, refrigerant heat exchanger 
(precooler, RHE) and mixture heat exchanger (solution heat exchanger, MHE). When 
these two heat exchangers were used in the system, the performance was increased. 
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Fig. 8. Exergy balance of a solar-driven ejector refrigeration system for the generating temperature of 90°C and 
evaporation temperature of 10°C [58]. 
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Fig. 9. Schematic of a absorption heat pump drawn by Sozen et al. [60]. 


For each individual component of the system shown in Fig. 9, the exergy destructions 
are written as follows [60]: 


EX¢est,cond = 77(EX7 — Exs) — Qoonall — (To/Tcona)], (52) 


Exdest,evap = m7(Ex19 — Ex11) + Devel > (To/Tevap)]. (53) 
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ExXdestabsor = 17 EX12 + rig6EX6 — rm EX1 — Oapsorll — (To/Tabsor)]s (54) 
ExX¢estpump = 71 (EX, — Ex2) + Wpumpe; (55) 
Ežaest,g = 113 EX3 — riigEX4 — rity EX; + O,[1 — (To/T»)], (56) 
EX¢estmix,HE = 112(EX2 — Ex3) + m4(Ex4 — Exs), (57) 
Exdest.HE = ?t7(Exg — Ežo + Ex, — Ei). (58) 


The exergetic coefficient of performance of the system for cooling and heating are as 
follows, respectively [60]: 
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Sozen et al. [60] performed the performance tests of the system under the climate 
condition of Ankara in Turkey. A parabolic collector was installed to obtain required 
temperatures. In the experiments, high temperature water obtained from the collector was 
used as heat source needed for the generator. The system design configuration was 
analyzed by using the experimental data. The effect of irreversibilities in thermal process 
on the system performance in the system were determined. Thermodynamic analysis 
showed that both losses and irreversibility had an impact on absorption system 
performance, while potential for the improvement in the components of the system was 
determined. The components of the evaporator and absorber of the system had a higher 
exergy loss than the other components. In order to improve the performance of the system 
thermally, those components should be modified. Typical COP. x, values of the system 
varied from 0.13 to 0.45. The highest and lowest COP.x, values were obtained at higher 
evaporator temperature (10°C) and the generator temperature of 75°C, and at 
Tevap = 10°C and Teona = 28 °C, respectively. 

Izquierdo et al. [61] studied on the entropy generated, the exergy destroyed and the 
exergetic efficiency of lithhum-bromide absorption thermal compressors of single and 
double effect, driven by the heat supplied by a field of solar thermal collectors. A schematic 
of a single effect thermal compressor fed by a thermal collector is illustrated in Fig. 10 [61], 
where the solar cases heat storage system was not included. It consisted of an absorber, a 
desorber, the heat recoverer, the solution pump and the throttling valve. The compressor, 
driven by the heat produced by a field of solar collectors, absorbed the low pressure. 

Izquierdo et al. [61] reported as follows: Two different applications, namely air-cooled 
and water-cooled units, were taken into account and compared with each other. Water- 
cooled compressors operated with temperatures and pressures lower than air-cooled 
compressors considering, in both cases, the same suction temperature, equal to 5°C. While 
the absorption temperature in water-cooled compressors could reach 40 °C, in air-cooled 
systems it could vary between 30°C and more than 50°C. Under these conditions, 
the discharge temperature (boiling temperature within the desorber) of a single effect 
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Fig. 10. Schematic of a single effect thermal compressor fed by a thermal collector drawn by Izquierdo et al. [61]. 


air-cooled unit ranged from 65 to 110°C, while the maximum discharge pressure was 
around 0.12 bar. The discharge temperatures (boiling temperature within the high pressure 
desorber) of the double effect air-cooled thermal compressor varied from 110°C for a final 
absorption temperature of 30-180°C for a final absorption temperature of 50°C. 
Discharge pressures could reach values of 0.3 and 1.5bar, respectively. The lithium- 
bromide air-cooled thermal compressors of double-effect could be viable with absorption 
temperatures around 50°C, when the temperature difference between the lithium-bromide 
solution and the outside air was about 8°C. The exergetic efficiency of the air-cooled 
compressor of double effect was found to be 10% greater than the exergetic efficiency of 
the air-cooled thermal compressor of single effect. The exergetic efficiency of the air-cooled 
thermal compressor of double effect was also determined to be 10% higher than the 
exergetic efficiency of the water-cooled thermal compressor of single effect. The exergy of 
the solar radiation was better used by the air-cooled thermal compressor as the 
temperature of the external fluid feeding the high desorber was higher. 

4.1.1.3.3. Adsorption units. Adsorption refrigeration is considered alternative to the 
conventional vapor compression refrigerator, especially in remote areas of the world 
without grid connected electricity. Activated carbon/methanol, activated carbon/ammo- 
nia, zeolite/water and silica gel/water are the adsorbent/adsorbate pairs commonly used in 
practically realized adsorption refrigeration cycles. A simple adsorption cycle indicating all 
the heat transfer for a complete cycle is shown in Fig. 11 [62]. Its principle of operation 
involves the processes of isosteric heating, desorption, isosteric cooling and resorption [62]. 

The thermodynamic efficiency (exergetic efficiency) relates the actual performance of a 
system to its expected maximum performance as determined from its ideal Carnot cycle. It 
may be written from a relation proposed by Pons et al. [63] as follows: 


COP, 
(Pa ma 
with 
l1 — Ta sor T 
COP max = ( ad / g) (62) 


(Tadsor/ Tevap) -1 
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Fig. 11. Schematic of a simple adsorption refrigerator [62]. 


and 


COP, = evap (63) 


where Ding and Os are the heat transferred during refrigeration and the heat used to 
generate refrigerant during generation, respectively [62]. 

Anyanwu and Ogueke [62] reported as follows: The COP and thermodynamic efficiency 
depended on the mass of adsorbate generated and the generation temperature, T,. Above a 
certain generator temperature, the mass of adsorbate generated was no longer a direct 
function of the temperature. In fact, in some cases, continued increase in generator 
temperature did not yield any adsorbate. Thus, an improved COP and consequently the 
thermodynamic efficiency would only come from an adsorbent that would readily desorb 
the adsorbate in its pores up to a low level of about 20—-25g of adsorbate per kg of 
adsorbent for the range of temperature obtainable using a conventional flat plate solar 
collector. This, however, was with the understanding that the refrigerant would have a 
high latent heat of vaporization. The activated carbon/ammonia for 0, —10 and —20°C 
evaporator temperatures and 47 °C condenser temperature gave thermodynamic efficiency 
values above 100% at generator temperatures of about 70°C. When the interest was air 
conditioning, the zeolite/water pair would be preferred while activated carbon/ammonia 
would be the preferred pair for other applications. 


4.1.1.4. Solar cookers. Solar energy, which is an abundant, clean and safe source of 
energy, is an attractive to substitute for the conventional fuels for cooking [33]. In the 
literature, there are two studies conducted on the exergetic evaluation of the performance 
of solar cookers. A simply designed and the low-cost parabolic-type solar cooker (SPC) 
shown in Fig. 12 was made and tested by Ozturk [33]. He experimentally assessed the 
energy end exergy efficiencies of the cooker. Petela’s analysis [32] was inspired by 
publication of Ozturk [33], in which for the first time the exergy efficiency was determined 
experimentally for the SPC of the cylindrical trough shape [32]. 


20 mm 
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wee profile) 


Fig. 12. A sketch of a solar parabolic cooker drawn by Ozturk [33]. 


The instantaneous exergy efficiency of the solar cooker can be defined as the ratio of the 
increased water exergy to the exergy of the solar radiation. Using Eqs. (21), (22) and (23b), 
it may be calculated as follows [33]: 


My Cy [(T wout = Twin) — To (In(Tw,out/ Tw,in)) | /At 
IpAcook [1 + 1/3(To/Ts)* — (4/3XTo/T9] ` 
Exergy efficiency of the SPC was also defined by Petela [32] as the ratio of the exergy of 
the useful heat Q;,,, at temperature 73, and of the exergy of solar emission at temperature 
T;. Using Eqs. (21), (22) and (24), it was proposed as follows [32]: 
Os,(1 — (To/T3)) 
Iy[1 + 1/3(To/Ts)* — (4/3)(To/Ts)] 


with the useful heat transferred through the wall of the cooking pot given by 


(64) 


Escook,1 = 


(65) 


Escook,2 = 


03, = —43U3(T3 — Tw), (66) 
where T3 is the absolute temperature of the surface 3, A3 the outer surface area 3 of the 
cooking pot, U; the heat transfer coefficient (which takes into account the conductive heat 
transfer through the cooking pot wall and convective heat transfer from the inner cooking 
pot surface to the water) and Tw the absolute temperature of water in the cooking pot, as 
shown in Fig. 13 [32]. 

Fig. 12 illustrates the dimensions of a SPC, which was exergetically tested by Ozturk 
[33]. The overall dimensions of the SPC are 0.90 x 1 m”. The SPC was constructed of steel 
profile and Cr—Ni alloy sheet. Structural frame of the SPC, made of steel profile of 
dimensions 0.15 x 0.15 m° and Cr-Ni alloy sheet, was screwed on the frame. The thickness 
of Cr-Ni sheet was 0.0005 m. In the center of the concentrating reflector, a cooking pot of 
0.10m width and 0.05m in depth was welded. The cooking pot was made of galvanized 
steel sheet of 0.005 m thickness. All outer surfaces of the cooking pot were painted with 
mat black. There were no clouds so the reflector acted only on direct solar radiation. The 
reflector could be adjusted by means of the two screws on both sides of the cooking pot. 
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Fig. 13. A schematic of a solar parabolic cooker drawn by Petela [32]. 


The emissivity of the cooking pot and shield surfaces were 0.87 and 0.45, respectively 
[33]. Assumed were the following data: the time interval Aż = 600s, the constant 
specific heat for water Cy = 4.2kJ/kg K and the solar radiation temperature T, = 6000 K. 
The experimental time period was from 10:00 to 14:00 solar time. During this 
period, it was found that the daily average temperature of water in the SPC was 333K 
and the daily average difference between the temperature of water in the cooking pot 
and the ambient air temperature was 31.6K. The exergy output of the SPC varied 
between 2.9 and 6.6 W. The exergy efficiency of the SPC was obtained to range from 0.4% 
to 1.25% [33]. 

Petela [32] derived equations for heat transfer between the three surfaces: cooking pot, 
reflector and imagined surface making up the system. His model allowed for theoretical 
estimation of the energy and exergy losses: unabsorbed insolation, convective and radiative 
heat transfer to the ambient, and additionally, for the exergy losses: the radiative 
irreversibilities on the surfaces, and the irreversibility of the useful heat transferred to the 
water. The exergy efficiency of the SPC was found to be relatively very low (about 1%), 
and to be about 10 times smaller than the respective energy efficiency. The reasons for this 
were mainly due to the large exergy of the escaping insolation and additionally due to the 
degradation of the insolation absorbed on the surfaces of the reflector and the cooking pot. 
The theoretically calculated values of the efficiencies were relatively close to the values 
experimentally measured by Ozturk [33]. The influence of the input parameters 
(geometrical configuration, emissivities of the surfaces, heat transfer coefficients and 
temperatures of water and ambience) was also determined on the output parameters, the 
distribution of the energy and exergy losses and the respective efficiencies [32]. 


4.1.1.5. Industrial process heat. Cylindrical parabolic collector systems look very 
promising for delivering industrial process heat applications in 95-350°C delivery 
temperature range [41]. In this regard, Eskin [42] made a simulation study on a solar 
process heat system shown in Fig. 14 [41], validated by experimental results and 
represented the unsteady performance analysis of this system. She also evaluated exergetic 
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Fig. 14. A schematic view of the solar process heating system [41] (I: cylindrical parabolic collector array, II: 
storage tank, III: circulating pump, IV: auxiliary heater). 


performance of the particular systems, while her analysis considered the unsteady-state 
thermal and exergetic assessment of the collector and the storage tank individually. 

The exergy efficiency of the cylindrical parabolic solar collector was calculated as 
follows: 


i 
Epar.cot = 1 — Ee (67) 
EXpar,s 
with 
EXpars = AW srad,be,d' 7 


The auxiliary heater in the system eliminated the difference between energy demand and 
solar energy harvesting. When the temperature of the water entering the load fell below the 
predetermined temperature level, the auxiliary heater was activated. The simulation results 
indicated that the exergetic efficiency (at a maximal value of about 37%) was highly 
dependent on the ratio of mass flow rates and the use of an auxiliary heater in the system. 


4.1.1.6. Solar desalination systems. There are very limited studies on the exergetic 
evaluation of solar desalination systems in the open literature. The Spanish research 
institution, Centro de Investigaciones Energtticas, Medioambientales y Tecnologicas and 
the German Deutsche Forschungsanstalt fur Luft- und Raumfahrt e.V. decided in 1987 to 
develop a solar thermal desalination project at the Plataforma Solar de Almeria. In July, 
1988, the first start-up of a solar multi-effect distillation system took place at the 
Plataforma Solar de Aimeria, a solar research centre located in southeastern Spain, near 
Aimeria. Its evaluation was finished in December, 1990. The plant, known as “Sol-14” 
was built and connected to the previous existing solar facilities as a result of a 
Spanish-German project [64]. The performance of this solar desalination system was 
evaluated by Garcia-Rodriguez and Gomez-Camacho [64] using energy and exergy 
analysis method in order to propose possible improvements of the system, as presented 
partly below. 

In the exergetic analysis of the system indicated in Fig. 15 [64], Eq. (15) was used and the 
results obtained on a component basis are listed in Table 3 [64]. It is obvious that solar 
multi-effect distillation with 14 cells had the lowest exergetic efficiency, while the thermal 
storage had the highest one. 


A. Hepbasli / Renewable and Sustainable Energy Reviews 12 (2008) 593-661 623 


Solar 
collector 
field 


12 15 13 


Fig. 15. Solar multi-effect distillation (MED) system at Phase I of the solar thermal desalination project [64]. (1) 
Solar input, (2) optic and thermal losses, (3—6) synthetic oils at various temperatures, (7) thermal losses, (8) 
saturated steam, (9) saturated water, (10) fresh water (product), (11) thermal losses, (12) seawater, (13) cooling 
outlet stream, (14) blowdown and (15) auxiliary power consumption. 


Table 3 

Exergetic performance of the equipment at Phase I of the solar thermal desalination project [64] 

Equipment F (kW) P (kW) È — P (kW) e= P/F (%) 
Solar collector 336.9 58.63 278.3 17.4 
Thermal storage 58.63 ~58.63 =0 ~ 100 

Boiler 58.63 26.21 32.42 44.7 
Solar multi-effect distillation (MED) 34.91 4.992 29.92 14.3 
Global system 345.6 4.992 340.6 1.4 


4.1.1.7. Solar thermal power plants. The solar thermal power system in general can be 
considered as consisting of two subsystems, namely, the collector—receiver subsystem and 
heat engine subsystem as illustrated in Fig. 16 [12]. The collector—receiver circuit consists 
of a number of parabolic trough collectors arranged in modules operating in tracking 
mode. The main components of the heat engine circuit are a heater (here the boiler heat 
exchanger), a turbine having two stages, a condenser, pump and a regenerator. Exergy 
analysis is conducted in two main subsystems along with their components, as given by 
Singh et al. [12]. 

The overall exergy efficiency of the collector—receiver circuit is calculated in a similar 
manner given by Eq. (22) as follows [12]: 


EXxucol- 
€0,col-rec = ees (69) 
Xpar,s 
with the useful exergy of the whole parabolic collector 
EXu col-rec = Npar,colEXu (70) 
and the exergy received by the parabolic collector given as: 
: . T 
Ežpars = Ô, (: = D] (71) 
S 


where Mpar,coi Is the number of the parabolic collectors, Ò, is the transferred solar heat and 
T; is the solar temperature. The exergy received by the parabolic collector may also be 


624 A. Hepbasli / Renewable and Sustainable Energy Reviews 12 (2008) 593-661 


Jp 
f, 
4 
[u] - 
pA 2 i 
Tfi t B 4 
' a 
17 VII 
C 5 
6 A 
10 
Solar collector receiver 


Lt- ý . | 
circuit - >- O | te 
8 9 


Heat engine circuit — pl 


Fig. 16. Schematic of a solar power thermal plant drawn by Singh et al. [12]. (I) Parabolic through collector, (II) 
boiler heat exchanger, (III) high-pressure turbine, (IV) low-pressure turbine, (V) condenser, (VI) valve, (VI) 
regenerator. 


calculated using Eqs. (25)-(26), and the relations given in Table 1. The useful exergy of the 
one parabolic collector (Ex,) may be found using Eq. (23b) or (23c). 

The total efficiency of the heat engine subsystem consisting of the boiler heat exchanger 
and the heat engine cycle is calculated as follows: 


W net 
= 2 
Eh,eng EXt eng (72) 
with the available exergy of the working fluid of the heat engine cycle given as: 
. . T 
Exheng = 0, (1 == £ j (73) 
Ran,0 


where Ò, is the useful transferred heat, Wnet is the net work done by the heat engine cycle 
and TRan.o is the ambient temperature of the Rankine cycle [12]. 

The overall exergy efficiency of the solar thermal power system is the product of the two 
separate efficiencies of the subsystems, as given below [12]: 


W vet 
EXpars 


Eh,eng = (74) 
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You and Hu [65] investigated the optimal thermal and exergetic efficiencies of a 
medium-temperature solar thermal power system using energy and exergy analysis 
method. A simple reheat-regenerative arrangement is schematically shown in Fig. 17 [65]. 

The exergetic efficiency of the cycle (including the exergy loss in the heating process) is as 
follows [65]: 


W turbine 


Exo (75) 


Ecycle = 
Exg is calculated using Eq. (23c) where Ty our and Ty,our are the temperatures of the 
solar heat carrier entering and exiting the heat exchanger(s), respectively. 
The exergetic efficiency of the collector is calculated from [65]: 


Exg 
pars = Ek, (16) 
with 
4 Tam 
Ex, = I7|1 -3T (1 — 0.28f) |. (17) 


The overall exergy efficiency of the combined system of the collector and the power 
cycle is the product of the two separate efficiencies of the subsystems, as given by the 
equation [65]. 

Esys = Ecycle€par,s- (78) 

Singh et al. [12] calculated the irreversibilities of the solar thermal power plant illustrated 
in Fig. 16 [12]. It was found that the collector and receiver had the highest irreversibilities 
and the efforts to reduce exergy loss should be made in this assembly. The component 
exergy efficiency values were as follows: Collector: 29.03%, receiver: 68.38%, collector- 
receiver: 19.85%, boiler heat exchanger: 96.42%, heat engine: 66.57% and overall system: 
12.74%. 

You and Hu [65] reported based their analysis that the reheat-regenerative arrangement 
was suitable for the medium-temperature solar thermal power generation purpose. There 
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Fig. 17. Schematic of one-stage reheat-regenerative power cycle and the parabolic trough collector drawn by You 
and Hu [65]. 
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was an optimum saturation temperature in the boiler or an optimum temperature the fluid 
entering the solar field (with the same exiting temperature) at which both the thermal 
efficiency and the exergetic efficiency of the combined system of the reheat-regenerative 
cycle and the solar collector hit the maximum values simultaneously. 


4.1.2. Photovoltaics 

Exergetic evaluation of PVs has been performed by some investigators [66,67] in hybrid 
(PV/thermal) systems as a part of the system. This evaluation is explained in Section 4.1.3 
in more detail. 

Electrical energy is not affected by ambient conditions and therefore is equivalent in 
work. If global irradiance is /, energetic efficiency of the solar cell is neg the instantaneous 
electrical exergy is then as follows [66]. 


Exe = Nscell = Escell!, (79) 


where Esce is the exergetic efficiency of the solar cell. 


4.1.3. Hybrid (PV/thermal) solar collectors 
The exergy efficiency of a hybrid solar collector, that generates both electric power and 
heat, may be calculated as follows [67]. 


Heonver! + ExXsolar 


€pv/Thermal = Ex (80) 
with 
s T uid — To - 
Exsolar = a Qsotar (81) 
fluid 
and 
Ex = 0.951, (82) 


where MNeonver is the conversion efficiency, J the global irradiance (W/m7), EXxsolar the exergy 
of heat (W/m), Ex; the exergy of global irradiance (W/m), Qj, the collected solar heat 
amount per unit time per panel area (W/m?) and Tayia the supply temperature of collector 
fluid (K). 

Fujisawa and Tani [66] defined the synthetic exergy of the PV/T collector as the total 
value of the electrical and thermal exergies as follows: 


Expy/t = Exe + Exthermal = (Ee F Ethermal) = épy/ TL (83) 
with 
Ei, = el (84) 
and 
; : -/Tt-T 
EXthermal = Oncarnot = o=) = Ethermall. (85) 


Fujisawa and Tani [66] designed and constructed a PV-thermal hybrid collector on their 
campus. The collector consisted of a liquid heating fiat-plate solar collector with mono-Si PV 
cells on substrate of non-selective aluminum absorber plate, The collector area was 
1.3 x 0.5m?. From the annual experimental evaluation based on exergy, Fujisawa and Tani 
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concluded that the PV/T collector could produce higher output density than a unit PV module 
or liquid heating fiat-plate solar collector. Using exergetic evaluation, the best performance of 
available energy was found to be that of the coverless PV/T collector at 80.8kWh/yr, the 
second to be the PV module at 72.6 kWh/yr, the third to be the single-covered PV/T collector 
at 71.5kWh/yr, and the worst to be the flat-plate collector at 6.0 kW h/yr. 

Saitoh et al. [67] performed field measurements from November 1998 to October 1999 at 
a low energy house at Hokkaido University. A system diagram of the experimental 
equipment is illustrated in Fig. 18 [67]. The power generated was measured by giving the 
variable load for the maximum power point tracking. Two hybrid solar panels were 
connected in parallel, while the brine was supplied at constant temperature by a fluid 
supply system with a circulating pump. The volumetric flow rate per panel was fixed at 1 1/ 
min. The tilt angle of the panel was 30°, which gave the maximum annual global 
irradiance. Except for winter, the mean conversion efficiency of array and the collector 
efficiency were stable at 8-9% and 25-28%, respectively. The dependency on solar energy 
was 46.3%. Energy and exergy efficiency values of single-junction crystalline silicon PV are 
illustrated in Table 4 [67]. The exergy efficiencies of solar collector, PV and hybrid solar 
collector were found to be 4.4%, 11.2% and 13.3%, respectively. Solar collector had the 
lowest value. In terms of energy quality, the hybrid solar collector was very effective [67]. 


4.2. Wind energy systems 


The exergy of fluid currents is entirely due to the kinetic energy of the fluid. The 
maximum work is extracted from a moving fluid when the velocity is brought to zero 
relative to the reference state. The kinetic energy of the fluid and its exergy have the same 
numerical value [14]. 


1 
Y fuid = 5P Vo (86) 


where p is the density of the wind and Vo is the wind speed. 
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Fig. 18. System diagram of experimental equipment drawn by Saitoh et al. [67]. 
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Table 4 
Comparison of energy and exergy efficiency values of solar collector, photovoltaic and hybrid collector [67] 


Type of efficiency Solar collector Photovoltaic Hybrid collector 


Energy efficiency (n, %) 


Heat 46.2 — 32.0 
Power — 10.7 10.6 
Total 46.2 10.7 42.6 
Exergy efficiency (e, %) 

Heat 4.4 — 2.1 
Power — 11.2 11.2 
Total 4.4 11.2 13.3 


Koroneos et al. [11] presented two diagrams where utilization of the wind’s potential in 
relation to the wind speed and exergy losses in the different components of a wind turbine 
(i.e., rotor, gearbox and generator) were included. They concluded as follows: “According 
to Betz’s law, wind turbine can take advantage of up to 60% of the power of the wind. 
Nevertheless, in practice, their efficiency is about 40% for quite high wind speeds. The rest 
of the energy density of the wind not obtainable is exergy loss. This exergy loss appears 
mainly as heat. It is attributed to the friction between the rotor shaft and the bearings, the 
heat that the cooling fluid abducts from the gearbox, the heat that the cooling fluid of the 
generator abducts from it and the thyristors, which assist in smooth starting of the turbine 
and which lose 1-2% of the energy that passes through them.” 

As for as studies conducted on exergetic assessment of wind energy systems, Sahin et al. 
[68] developed a new exergy formulation for wind energy, which was more realistic and 
also accounted for the thermodynamic quantities enthalpy and entropy. The relation 
developed for the total exergy for wind energy is as follows, while the input and output 
variables for the system considered are illustrated in Fig. 19 [68]. 


T Py È 
Ex = Egen + CT — T1) + mT at (<, ia = R Iin — 9) (87) 
Tı Pi Tu 


with 
Oross = mCp(Tat = Taye), (88) 


where Ea is the generated electricity (i.e., the total kinetic energy difference), T,, is the 
atmospheric temperature, P is the pressure at the exit of the wind turbine for a wind speed 
V, and P; is the pressure at the inlet of the wind turbine for a wind speed Vi, Q,,., 
represents heat losses of wind turbine and Tave is the mean value of input and output wind 
chill temperatures. 

Sahin et al. [68] reported as follows: “The exergy analysis of wind energy shows that 
there are significant differences between energy and exergy analysis results. According to 
one classical wind energy efficiency analysis technique, which examines capacity factor, the 
resultant wind energy efficiency is overestimated relative to what it really is. The capacity 
factor normally refers to the percentage of nominal power that the wind turbine generates. 
Average differences between energy and exergy efficiencies are approximately 40% at low 
wind speeds and up to approximately 55% at high wind speeds”. 
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Fig. 19. Wind turbine and representative wind energy input and output variables drawn by Sahin et al. [68]. 


Ozgener and Ozgener [69] performed exergy analysis of a wind turbine system (1.5 kW) 
located in Solar Energy Institute of Ege University (latitude 38.24N, longitude 27.50E), 
Izmir, Turkey. They reported that exergy efficiency changed between 0% and 48.7% at 
different wind speeds based on a dead state temperature of 25°C and a atmospheric 
pressure of 101.325 kPa, considering pressure differences between state points. Depending 
on temperature differences between state points, exergy efficiencies were found to be 
0-89%. 

The exergetic efficiency of a turbine is defined as a measure of how well the stream 
exergy of the fluid is converted into useful turbine work output or inverter work output. By 
using inverter work output as useful work and applying this to wind turbine, the following 
relations are obtained [69]. 


We io 
p = — 9 
E= Ex, — E% m 
with 
EXdest = (EX1 — Exo) — We, (90) 


where W cio is the power at inverter output and Ex¢es is the exergy destruction, while Ex; 
and Ex, are exergy rates calculated using Eqs. (12) and (13). 


4.3. Geothermal energy systems 
The specific exergy of a geothermal site depends on the temperature at the collector, the 


pressure and composition of fluid in the reservoir, and the ability of the surrounding rock 
or fluid to transport thermal energy to the collector. Utilization of a geothermal site often 
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changes the thermal or fluid properties of the reservoir to a new equilibrium dependent on 
the method of energy extraction [14]. 


4.3.1. Classification of geothermal resources by exergy 

Geothermal resources have been classified as low, intermediate or high enthalpy 
resources according to their reservoir temperatures. The temperature ranges used for these 
classifications are arbitrary and they are not generally agreed upon. Temperature is used as 
the classification parameter because it is the earliest to measure and understand. In 
addition, temperature or enthalpy alone can be ambiguous defining a geothermal resource 
because two independent thermodynamic properties are required to define the thermo- 
dynamic state of fluid. Geothermal energy is already in the form of heat, and from the 
thermodynamic point of view, work is more useful than heat because not all the heat can 
be converted to work. Therefore, geothermal resources should be classified by their exergy, 
a measure of their ability to do work [70,71]. 

Lee [71] proposed a new parameter, the so-called specific exergy index (SExI) for better 
classification and evaluation as follows: 

Morine — 273.1 65prine 


SExI = aCe ; (91) 


which is a straight line on an h—s plot of the Mollier diagram. Straight lines of SExI = 0.5 
and SExI = 0.05 can therefore be drawn in this diagram and used as a map for classifying 
geothermal resources by taking into account the following criteria: 


è SExI<0.05 for low-quality geothermal resources, 
è 0.05<SExI<0.5 for medium-quality geothermal resources, and 
è SExI>0.5 for high-quality geothermal resources. 


Here, the demarcation limits for these indices are exergies of saturated water and dry 
saturated steam at 1 bar absolute. 

In order to map any geothermal field on the Mollier diagram as well as to determine the 
energy and exergy values of the geothermal brine, the average values for the enthalpy and 
entropy are then calculated from the following equations [72]: 


Morine = ee, (92) 
Sonae = Spel, (93) 

The above classification has been used by some investigators in their studies. In this 
regard, Quijano [72] performed exergy analysis of the Ahuachapan and Berlin geothermal 
fields. Plotting the thermodynamic conditions of both fields on the Molliere diagram and 
computing the specific exergy index (SEI), the fields were classified on the medium high 
exergy zone. Ozgener et al. [73] and Baba et al. [74] applied the SExI to the Balcova 
geothermal field, which is located in the western part of Turkey, and is endowed with 
considerably rich geothermal resources. The number of the wells in operation in this field 
may vary depending on the heating days and operating strategy. Taking into account the 
eight productive wells, which were in operation at the date considered and using Eq. (91) 
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together with the values given in Table 5 [74], the SExI was found to be 0.07. This 
represented that the field fell into the medium quality geothermal resources according to 
the classification of Lee [71]. For classifying the Balcova geothermal field according to its 
reservoir, the weighted average value of the temperature (Thrine) was obtained to be 
113.9 °C by using the similar relation given in Eq. (91). This field fell into the intermediate 
enthalpy resources according to the classifications of Muffler and Cataldi [75] and 
Benderitter and Cormy [76], while it is characterized as low enthalpy resources according 
to the classification of Hochstein [77]. This clearly indicated that the classification with 
reference to SExI was more meaningful as there was no general agreement on the arbitrary 
temperature ranges used in the classification of geothermal resources by temperature. 
Ozgener et al. [78,79] also determined the SExI values for the Gonen and Salihli 
geothermal fields in Turkey as 0.025 and 0.049, respectively. They concluded that both 
these fields fell into the low quality geothermal resources according to the classification of 
Lee. Etemoglu and Can [80] classified geothermal resources in Turkey based on SExI. In 
their analysis, the triple point of water was selected as the dead-state conditions since the 
enthalpy and entropy values were zero. SExI values of geothermal resources over 0.56 were 
classified as high exergy resources, SER <0.05 low exergy and between of these values were 
medium exergy resources. 


4.3.2. Direct utilization 

4.3.2.1. District heating. Although various studies have been undertaken to perform 
energy and exergy analyses of geothermal systems, the exergy analysis of geothermal 
district heating systems (GDHSs) based on the actual data was performed by Ozgener 
et al. [81] for the first time (to the best of authors’ knowledge). In this context, the first 
study included a comprehensive energy and exergy analysis conducted on the Balcova 
geothermal district heating system located in the town Balcova, Izmir, Turkey, using actual 
data sets obtained as of January 1, 2003 for system analysis, performance evaluation and 
efficiency improvement. Since then, Ozgener et al. have conducted various studies on 
exergetic [82,83] and exergoeconomic analyses [84,85] of GDHSs, namely Balcova, Salihli 
and Gonen GDHSs installed in the provinces of Izmir, Manisa and Balikesir, Turkey 


Table 5 
Some thermodynamic properties of the wells in the Balcova geothermal field in Izmir, Turkey (as of 1 January 
2003) [74] 


Code of well Temperature T Mass flow rate 7 Pressure P Specific enthalpy Specific entropy 
CC) (kg/s) (kPa) h (kJ/kg) s (kJ/kg K) 

BD2 131.6 18.31 413.4 552.98 1.6507 
BD3 120 18.65 522.83 503.72 1.5270 
BD4 135 25.32 539.04 567.57 1.6863 
BD5 111.4 3.74 607.95 467.35 1.4333 
BD7 110.6 20.01 496.49 463.89 1.4245 

B4 98.8 13.16 368.81 413.99 1.2926 

B5 103.6 23.10 332.34 434.20 1.3467 
B10 96 25.90 364.76 402.19 1.2607 
BD8 66.1 112.02 474.20 276.8 0.9061 
(reinjection) 


Average values 113.9 148.19 — 477.87 1.4578 


632 A. Hepbasli / Renewable and Sustainable Energy Reviews 12 (2008) 593-661 


as well as on investigating the effect of dead state temperature on the energy and 
exergy efficiencies of GDHSs [86,87]. They also performed a key review on the 
performance improvement aspects of GDHSs [88] and prepared a general paper on 
geothermal energy including the exergetic evaluation of GDHSs for the Encyclopedia of 
Energy Engineering [89]. 

The following relations have been summarized from the previous studies performed by 
Ozgener et al. [81-91] as follows: 

A schematic diagram of the simplified GDHS, where user system is heated by 
geothermal energy, is illustrated in Fig. 20 [88]. As can be seen in this figure, a GDHS 
consists mainly of three circuits, namely (a) energy production circuit (EPC; geothermal 
well loop), (b) energy distribution circuit (EDC; district heating distribution network, and 
(c) energy consumption circuit (ECC). In evaluating the system performance from the 
exergetic point of view, the two items of the whole system play a big role. These include: (i) 
the geothermal brine exergy inputs from the production field and (ii) exergy destructions in 
the system taking place as the exergy of the fluid lost in the heat exchanger, the natural 
direct discharge of the system (pipeline losses), and the pumps. 

For a general steady-state, steady-flow process, the two balance equations, namely mass 
and exergy balance equations, may be written as follows: 

Using Eq. (1), the mass balance equations for two various GDHSs, namely GDHSs with 
reinjected thermal water (i.e., Balcova GDHS) and without reinjected thermal water 
(i.e., Salihli GDHS), may be written as follows, respectively [88]: 


n 
> My T — M, — mq = 0 (94) 
i=1 
and 
n 
> Mw T — Ma = 0, (95) 
j=l 
EPC Central Heating 
Station 
EDC 
a) Ener; roduction 
, cycle (EPC) (b) Energy distribution ECC 
cycle (EDC) > 
O (c) Building substations 
a z 
za] (Energy consumption 
EDC EDC circuit -ECC-) 
@) BCC ECC 


EPC Circulating pump Circulating pump 


H, 


Reinjection well 


Fig. 20. Basic simplified geothermal district heating energy system schema [90]. 


A. Hepbasli / Renewable and Sustainable Energy Reviews 12 (2008) 593-661 633 


where my 7 is the total mass flow rate at wellhead, m, is the flow rate of the reinjected 
thermal water and mqa is the mass flow rate of the natural direct discharge. 

The exergy destructions in the heat exchanger, pump and the system itself are calculated 
as follows [88]: 


ÈXdest, HE = EXxin = ÈXout = ÈXdest, (96) 

ÈXdest pump = W pimp = (EXout = Exin) (97) 
and 

EX dest;systei = 5 ÈXdest HE F 5 Exdestpump: (98) 


Using Eq. (15), the exergy efficiencies of two various GDHSs with reinjected thermal 
water and without reinjected thermal water are found from the following relations, 
respectively [88]. 

EXuseful HE 1 ËXdest,system + EXreinjected + ÈXnatural discharged 


EGDHS,wr = = ; (99) 
EXprine Exprine 


ExXuseful, HE =] EXdest,system F EXnatural discharged 


> (100) 


é = z z 
coven EXprine EXprine 
where Ex¢est system iS the exergy destructions in the system taking place as the exergy of the 
fluid lost in the heat exchanger and the pumps, EXreinjectea is the exergy loss due to the 
reinjected water and EXpatural discharged 18 the exergy destruction due to the natural direct 
discharge of the system (due to a significant amount of water leaks; pipeline losses). 

Table 6 [81—83,86,87,90,91] shows a comparison of exergy efficiency values for GDHSs 
installed in Turkey. As can be seen in this table, exergy efficiencies range from 42.94% to 
66.13% depending on the various dead state temperatures. As expected, the lower the 
reference state temperature, the significantly larger the exergy losses in the system pipeline 
and heat exchangers. However, exergy rates of pipeline losses of the system, heat 
exchangers, and useful exergy increase considerably. The reason for this rapid rise in 
exergy rates is due to a decrease in the ambient temperature [88]. 

Ozgener et al. [88] also studied on the effect of varying dead state temperatures on the 
exergy efficiencies of the whole GDHS. Exergy efficiency values were correlated with high 
correlation coefficients for the Balcova, Gonen and Salihli GDHSs. Of these, the 
correlation for the Balcova GDHS is given below [87] 


Bakova, GDHS = 0.000006 — 0.0000573,, — 0.1279 Tam + 45.21, (101) 


where Tam is the ambient temperature in °C. 


4.3.2.2. Ground-source (geothermal) heat pumps. Ground-source heat pumps (GSHPs) 
have been used for years in the developed countries due to their higher energy utilization 
efficiencies compared with conventional heating and cooling systems. There are two main 
types of GSHP system. (a) ground-coupled (vertical or horizontal)-closed loop and (b) 
water source-open loop. Ground-coupled heat pumps are known by a variety of names. 
These include GSHPs, earth-coupled heat pumps, earth energy heat pumping systems, 
earth energy systems, ground-source systems, geothermal heat pumps, closed-loop heat 


Table 6 
Comparison of various geothermal district heating systems in terms of exergy destructions and exergy efficiencies [8 1—83,90,91,86,87] 
Name of Location/ Dead state Date of data Total exergy Exergy destructions/losses (% of the total (brine) exergy Exergy References 
GDHS country temperature used input (kW) input) efficiency 
CC) (%) 
Pumps Heat Thermal Natural 
exchangers reinjection direct 
discharge 
(pipe line) 
Balcova Izmir/Turkey 13.1 1 January 9164.29 3.22 7.24 22.66 24.1 42.94 81] 
GDHS 2003 
Gonen Balikesir/ 6 1 February 2657.5 14.81 711 12.96 1.06 64.06 82] 
GDHS Turkey 2004 
Salihli Manisa/ 2.9 1 February 2564 2.22 17.90 — 20.44 59.44 83] 
GDHS Turkey 2004 
Balcova Izmir/Turkey 13.1° 2 January 14808.15 1.64 8.57 14.84 28.96 46.00 84] 
GDHS 2004 
Gonen Balikesir/ 10° 1 February 2333.33 17.45 8.19 7.37 0.86 66.13 85] 
GDHS Turkey 2004 
Salihli Manisa/ 25° 1 February 1524 4.27 31.17 = 8.98 55.58 88] 
GDHS Turkey 2004 
Balcova Izmir/Turkey 11.4% 2 January 14390 1.74 8.83 14.18 28.7 46.55 89] 
GDHS 2004 


“The study was conducted for four various dates , namely 1 January 2003 , 2 December 2003, 2 January 2004 and 2 February 2004 at dead state temperatures of 


13.1, 12, 10.4 and 6.5 °C, respectively. 


>A parametric study was performed o investigate how varying reference how varying reference temperature from 0 to 20°C would affect the exergy efficiencies. 

“The study was conducted for 1 December 2003, 1 January 2004 and 1 February 2004 at dead state temperatures of 10.8, 9.7 and 2.9, respectively. In addition, two 
parametric expressions of exergy efficiencies were developed as a function of the reference temperature, by varying the dead state temperature from 0 to 25°C. 

“The reference-state temperature value of 11.4°C was obtained from the average ambient temperature data measured between 2000 and 2004 for second January in 
Izmir. The dead state temperature was considered varying from 0 to 25°C to conduct a parametric study. 


veo 
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pumps, solar energy heat pumps, geoexchange systems, geosource heat pumps, and a few 
other variations [92]. 

Fig. 21 [93] illustrates a schematic diagram of a GSHP system in the heating mode as an 
illustrative example. The system was designed and installed in the Solar Energy Institute of 
Ege University, Izmir, Turkey. The system was commissioned in May 2000, while GSHP 
systems have been applied to the Turkish residential buildings since 1997. 

The performance tests have also been performed since then. This system has been 
recently integrated with a solar collector and a greenhouse [94]. It mainly consists of three 
separate circuits: (i) the ground coupling circuit (brine circuit or water-antifreeze solution 
circuit), (ii) the refrigerant circuit (or a reversible vapor compression cycle) and (iii) the fan 
coil circuit (water circuit). 

In Table 7 [95], Eqs. (102)-(126) used in modeling and analyzing GSHP systems are 
summarized [95]. These relations have been obtained for each of the GSHP components 
illustrated in Fig. 21 using mass, energy, entropy and exergy balance equations as well as 
the exergy destructions obtained using exergy balance equations given in Section 3). 


Il 


7 ST T 8 
© Ground level 


VI 


Fig. 21. The main components and schematic of a ground-source heat pump system [93]. (1) Compressor, 
(ID) condenser, (HI) capillary tube, (IV) evaporator, (V) fan-coil unit, (VI) ground heat exchanger. 


Table 7 


Mass, energy, exergy and entropy balance relations of a GSHP system [95] 


Toit 


No Element Mass analysis using Eq. (1) (the The irreversibility or the exergy destroyed using Eq. The irreversibility or the exergy destroyed using 
conversion of mass principle) (5a) (the exergy balance) entropy balance and Eq. (10) (the decrease of exergy 
principle) 
I Compressor my = iy = Ty (102) Ex, + Weompe = xaa + hi 03) Íi = Torin (sı — 52a) (105) 
Ty = h(i — oq) + Weomp 04) 
Il Condenser m, = tin =m, (106) Exo, + Exg = Èx; + Exs + fy 08) İn = Toltine(s3 — 52a) + titw(ss — s6)] (110) 
Tis = m6 = Mw (107) İn = rite (hog — Wy) + tite (Wg — Ws) 09) 
Ill Throttling valve m3 = Mg = Mr (111) Ex; = Exy + İm 12) Tin = Tom (s4 — s3) (114) 
(capillary tube) 
Tin = ri. (Wy — Wy) 13) 
IV Evaporator my = m = iy (115) Exy + Exg = Ex; + Ex7 + İiv 16) hy = Tolin (sı — s4) + titwa(s7 — 5g)] (118) 
ly = Tie (Wg — Wy) + MitwalWg — Wr) 17) 
v Fan coil units Mairin = Mairout = Mair (119) Exs = Exo + On (=s = T) + iy 20) iv=T |i ecija Qr l (122) 
in,air Tinair 
A : i Ti 
iy = tls We) ~ Oe(1 = 4 ) 21) 
in,air 
VI Ground heat m7 = Mg = Mya (123) r- ; Toit — To _f j š ; 0. 
exchanger Ex7 + Qon To Paek i 9 Tvr = Po) mwa(ss — 87) — Tai (126) 
A : ; 7 
ha = tials =) + Oah 1-7) 25) 


99 
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For further details, some relevant studies published should be taken into account 
[93-97]. 

Using Eq. (15), exergy efficiency relations for the GSHP unit and the whole system are 
written as follows: 


EXxheat D EXcond,in = ExXcond,out 


&0,R,HP = = - : (127) 
° W ain W ain 
EXheat 
eae - - l (128) 
cael W ain T W pumps F W fan-coil 
COP, h 
= 129 
Ene = SOPor” (129) 
Þr HP 
EHP = =~, (130) 
Fryp 
P 
go. , (131) 
FT sys 


Exergy analysis has been applied to various types of heat pumps by a number of 
researchers [i.e., 41,50,51,56,98—-100]. However, the studies on exergy analysis of GSHPs are 
relatively few in numbers [93—97,101—103]. Some of them [i.e., 43,104] have been included in 
the exergoeconomic (a combination of exergy and economics) analysis of GSHP systems. The 
first detailed exergetic evaluation of a GSHP system was made by Hepbasli and Akdemir [95] 
to the current knowledge of the author, while Ozgener and Hepbasli [105] also conducted a 
comprehensive review on the energy and exergy analysis of solar assisted heat pump systems. 
In this regard, the performance of this GSHP system with a 50m vertical 32 mm diameter U- 
bend ground heat exchanger was assessed using energy and exergy analysis methods. Fig. 22 
[94] illustrates a graphical representation of the exergy balance, namely the Grassmann 
diagram, for the GSHP system. This diagram, on the other hand, gives quantitative 
information regarding the proportion of the exergy input to the system that is dissipated in 
each of the components. As can be seen in Fig. 22, the highest irreversibility occurred in sub- 
region I, the motor-compressor subassembly. The total magnitude of the losses was over 56% 
of the actual power input, while the mechanical-electrical losses accounted for 46% of that. 
These losses were due to the electrical, mechanical and isentropic efficiencies and emphasize 
the need for paying close attention to the selection of this type of equipment, since 
components of inferior performance could considerably reduce the overall performance of the 
system. The second largest irreversibility was due to the condenser. This was partly due to the 
large degree of superheat achieved at the end of the compression process, leading to large 
temperature differences associated with the initial phase of heat transfer. The third highest 
irreversibility was in the capillary tube due to the pressure drop of the refrigerant passing 
through it. Besides this, the evaporator had the lowest irreversibility on the basis of the heat 
pump cycle [95]. 

Table 8 [93-97,103,106,107] lists various ground source (geothermal) heat pumps 
analyzed using exergy analysis method Five types of exergy efficiencies given in 
Eqs. (127)-(131) were used in the calculations. It is obvious from this table that the 
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100% 


4.81% 


Fig. 22. Exergy diagram (the Grassmann diagram) for the GSHP system [94]. (I) Compressor, (II) condenser, 
(II) capillary tube, (IV) evaporator, (V) fan-coil units, (VI) ground heat exchanger. 


exergy efficiency values of GSHP units vary between about 67% and 92% on the product/ 
fuel basis depending on the various dead state temperatures. The values for the exergy 
improvement rate and thermodynamic parameters calculated using Eq. (16) and 
Eqs. (17)-(20) were also included in this table. 


4.3.2.3. Greenhouses. Good plant-growth conditions can be achieved by using green- 
houses. A greenhouse can be managed to protect the plants by creating a favorable 
environment, allowing intensive use of soil, and helping sanitary plant control. A variety of 
heating systems are being used in present day greenhouses to meet the heating and cooling 
requirements. As an alternative, there is great potential in employing a heat pump system 
for greenhouse air-conditioning based on its ability to perform the multi-function of 
heating, cooling and dehumidification [11]. 

Although various studies were undertaken to evaluate the performance of solar assisted 
heat pump systems, to the best of author’s knowledge, two studies [11,14] have appeared in 
the open literature on the experimental performance testing of a solar assisted vertical 
GSHP system for greenhouse heating. The first one [11] was related to an exergetic 
assessment at a fixed dead state temperature, while the second one [14] included a 
parametric study investigating the effect of varying dead state temperature on the exergy 
efficiencies. 


Table 8 
Various ground source (geothermal) heat pumps analyzed using exergy analysis method [93-97,103,106,107] 
Location Type of Capacity Dead Date Exergy destructions (rate of improvement potential) (% of the whole heat pump system) Eq. no. Exergy Thermodynamic Ref. 
of the heat pump/ in the state of data for efficiency parameters for 
system operating mode temp. used Ground source heat pump system (Whole system)? exergy (%) compressor (% 
built/ mode operated (°C) eff. of the whole system) 
country (kW) 
Ground source heat pump unit Other components Heat Overall 7 ô č f 
pump 
I Com. II Con. III Thr. IV V Fan-coil/ VI GHE VII Cir. Vill 
valve Evap. Con. WT /PHEx pumps Sol. 
col. 
Izmir/ Vertical 3.41 20 Feb. 2001 16.54 27.91 21.62 7.83 21.14 (127) 2.94 51.8 23.6 35.4 40.2 [95] 
Turkey GSHP/heating (129) 3.84 
(for space 
heating) 
Izmir/ Vertical 4.27 2.2 7 Jan.2004 51.83 8.00 8.36 9.42 17.67 (127) 9.07 [93] 
Turkey | GSHP/heating (57.12) (3.03) (2.43) (14.30) (17.21) 
(for space (128) 8.4 
heating) (130) 66.8 
(131) 66.6 
Izmir/ Solar assisted 3.98 10.93 17 Dec. 28.33 13.85 11.34 8.19 30.22 2.52 1.83 and 0.63 (130) 71.8 28.3 9.1 13.5 30.4 [94] 
Turkey vertical GSHP/ 2003-17 3.09 (131) 67.7 
heating March 2004 
(for greenhouse 
heating) 
Izmir/ Solar 4.15 —0.69 7 Jan. 2004 23.06 10.45 6.53 9.79 43.11 2.61 1.24 and 1.30 (130) 76.2 23.1 5.6 74 21.1 [96] 
Turkey assisted 1.89 (131) 75.6 
vertical GSHP/ 
heating 


(for greenhouse 
heating) 
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Table 8 (continued) 


Location Type of Capacity Dead Date Exergy destructions (rate of improvement potential) (% of the whole heat pump system) Eq. no. Exergy Thermodynamic Ref. 
of the heat pump/ in the state of data for efficiency parameters for 
system operating mode temp. used Ground source heat pump system (Whole system)" exergy (%) compressor (% 
built/ mode operated (°C) eff. of the whole system) 
country (kW) 
Ground source heat pump unit Other components Heat Overall 7 ô id f 
pump 
I Com. II Con. Ill Thr. IV V Fan-coil/ VI GHE VII Cir. VII 
valve Evap Con. WT /PHEx pumps Sol. 
col. 
Izmir/ Solar assisted 3.97 1.00 17 Dec. 3.13 17.70 13.68 6.43 49.08 8.61 0.8 0.57 (130) 91.80 [103] 
Turkey vertical GSHP/ 2003-17 (131) 86.13 
heating (for March 2004 
greenhouse 
heating) 
Erzurum/ GSHP with 7.02 25 45.92 15.83 6.27 11.69 14.64 /5.66 (130) 4.64 [97] 
Turkey water source/ (131) 1.44 


heating (for 
space heating)” 


“I-Com: compressor, I-Con: condenser, II-Thr. valve: throttle (expansion) valve, IV-Evap: evaporator, V: Fan-coil,/Con. WT: condenser water tank, VI-GHE: 
ground-heat exchanger/PHEx: plate heat exchanger, VII-Cir. Pumps: circulating pumps, VIII-Sol. col.: solar collector. 
>This system was designed and constructed by Kara [106] in 1999 on the base of a PhD thesis in the Mechanical Engineering Department, Ataturk University, Erzurum 
in Turkey in order to evaluate geothermal resources with low temperatures. The performance of the system was evaluated by Kara and Yuksel [107] using the energy 


analysis method only, while this system was evaluated by Akpinar and Hepbasli [97] using exergy analysis method. 
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A schematic diagram of the constructed experimental system is illustrated in Fig. 23 
[11,14]. This system mainly consists of three separate circuits as follows:(i) the ground 
coupling circuit with solar collector (brine circuit or water-antifreeze solution circuit), (ii) 
the refrigerant circuit (or a reversible vapor compression cycle) and (iii) the fan-coil circuit 
for greenhouse heating (water circuit). 

The exergetic results of these two studies are illustrated in Table 8. The highest 
irreversibility on a system basis occurred in the greenhouse fan-coil unit, followed by the 
compressor, condenser, expansion valve and evaporator, sub-regions I and V for the 
GSHP unit and the whole system, respectively. Besides this, the remaining system 
components have a relatively low influence on the overall efficiency of the whole system. 
Experimental have also showed that monovalent central heating operation (independent of 
any other heating system) could not meet overall heat loss of greenhouse if ambient 
temperature was very low. The bivalent operation (combined with other heating system) 
could be suggested as the best solution in Mediterranean and Aegean regions of Turkey, if 
peak load heating could be easily controlled [11]. 


4.3.2.4. Drying. Large quantities of food products are dried to improve shelf life, 
reduce packing costs, lower shipping weights, enhance appearance, encapsulate original 
flavor and maintain nutritional value. In evaluating the performance of food systems, 
energy analysis method has been widely used, while studies on exergy analysis, especially 
on exergetic assessment of drying process, are relatively few in numbers [108-112]. In 
these previous studies, the drying process was thermodynamically modeled by Dincer 
and Sahin [109], while drying of different products such as pistachio [108], red pepper 
slices [110], potato [111], pumpkin [112], wheat kernel [113] and apple slices [114] was 
evaluated in terms of energetic and exergetic aspects using various drying devices, such as a 
fluidized bed dryer, a solar drying cabinet, and cyclone type dryers [115]. Based on the 
previous studies, the first study on exergetic evaluation of a drying process using a GSHP 
drying system appears to be done by Kuzgunkaya and Hepbasli [115], who used laurel 
leaves as a product being dried, and made a comparison of exergetic efficiency values 
obtained from various studies on food drying systems. The main objective of this 
contribution was to perform exergy analyses of the drying process in terms of drying of 
laurel leaves under different operating conditions for the assessment of the drying 
performance. 

Fig. 24 [115] illustrates a schematic diagram of this system, which consists of mainly 
three separate circuits, namely: (I) the ground coupling circuit (brine circuit or 
water—antifreeze solution circuit), (II) the refrigerant circuit (or a reversible vapor 
compression cycle) and (II) the drying cabinet circuit (air circuit). The main components 
of the heat pump system are an evaporator, a condenser, a compressor and an expansion 
valve. To avoid freezing the water under the working condition and during the winter, a 
10% ethyl glycol mixture by weight was prepared. The refrigerant circuit was built on the 
closed loop copper tubing. The working fluid was R-22 [115]. 

There are mainly two ways of formulating exergetic efficiency for drying systems 
[108,113,115]. The first one can be defined as the ratio of the product exergy to exergy 
inflow as follows [108,115]: 


_ Exergy inflow — Exergy loss _ l Exioss (132) 
a Exergy inflow 7 EXin 
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Fig. 23. A schematic diagram of a solar asissted vertical ground source heat pump system for greenhouse heating 
[11,14]. 
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Fig. 24. Schematic diagram of the heat pump tray drying system [116]. (1) Refrigerant inlet to the expansion 
valve, (2) refrigerant outlet from the compressor, (3) water inlet to the ground heat exchanger, (4) water outlet 
from ground-heat exchanger, (5) air inlet to the condenser, (6) air outlet from the condenser, (7) air outlet from 
the drying chamber. 
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The second one may be defined on the product/fuel basis, as given in Eq. (133). The 
product is the rate of exergy evaporation (Exevap) and the fuel is the rate of exergy drying 
air entering the dryer chamber (Èxaa). In this regard, exergy efficiency may be written as 
follows [113]. 


ÈXevap 
= 133 
Nex = Ena (133) 
with 
é To \ - 
Pea (1 = =) Osas| (134) 
and 
O evap = My Neg, (135) 


where Oe is the heat transfer rate due to phase change, T,n, is the exit temperature of the 
material, 77, is the mass flow rate of water and Afg is the vaporization latent heat. 

The laurel leaves were sufficiently dried at the temperatures ranging from 40 to 50°C 
with relative humidities varying from 16% to 19% and a drying air velocity of 0.5 m/s 
during the drying period of 9h. The exergy efficiency values were obtained to range from 
81.35% to 87.48% using Eq. (132) based on the inflow, outflow and loss of exergy, and 
9.11% to 15.48% using Eq. (133) based on the product/fuel basis between the same drying 
air temperatures with a drying air mass flow rate of 0.12 kg/s [116]. 

A comparison of exergy efficiency values in drying various products is made in Table 9 
[108,112,113,115]. It is obvious from this table that the exergy efficiency values calculated 
using Eq. (132) are quite higher than those using Eq. (133). It may be concluded that it is 
more meaningful, objective and useful to assess the performance of the drying process 
relative to the performance of similar drying processes on the product (or benefit)/fuel 
basis [116]. 


4.3.3. Indirect utilization (geothermal power plants) 

Commercial delivery of geothermally generated electric power occurred in 1914 when a 
250kW unit at Larderello provided electricity to the nearby cities of Volterra and 
Pomarance [116]. Since then, many geothermal power plants have been installed. The 
performance of these systems has been mostly evaluated using energy analysis method, 
while the number of the geothermal plants, of which performance has been assessed using 
exergy analysis method, is very low. The concept of exergy was first used to analyze a 
geothermal power plant by Badvarsson and Eggers [117]. Their illustrative example, which 
compared the performances of single and double flash cycles based on a reservoir water 
temperature of 250°C and a sink condition of 40°C, gave exergetic efficiencies of 38.7% 
and 49%, respectively, assuming 65% mechanical efficiency [71]. 

Using Eq. (15) along with its concept and applying it to a power plant as a whole, the 
overall exergetic efficiency reduces to a very simple formula, namely, the ratio of the net 
power output (W net,plant) to the exergy of the motive fluid serving as the energy source for 
the plant [118]. That is: 

W net,plant , (136) 


Eg.pp = ExXwh 
W. 


Table 9 
Comparison of exergy efficiency values calculated by some investigators in drying various products [108,112,113,115] 
Investigators/year of | Type of dryer Product dried Air Exergy efficiency 
publication 
Type Initial Final Velocity (m/s) Temperature Relative Eq. no. Value (%) 
moisture moisture CC) humidity (%) 


content (%) content on 
dry weight 
basis (%) 


Syahrul et al. [113] Fluidized bed Red-spring 31.7-32.6 15 1.95 40.2-65 18.5-21.1 (133) 4-12 
dryer wheat 
Corn 24.6-25.6 2.22-2.24 50-63 15.2-17.5 2-16 
Midilli and Kucuk Solar drying Shelled/ 26.95/29 1.23 40-60 37-62 (132) 10.86-100 
[108] cabinet unshelled 
pistachios 
Akpinar et. al. [112] Cyclone type Pumpkin 6 1.5 60 10-20 (132) 30.81—100 
dryer slices 
6 70 37.89-100 
4 80 46.97-100 
Kuzgunkaya and Ground- Laurel leaves 94.4 12 0.5 40 16-19 (132) 81.35-87.48 
Hepbasli [115] source heat 50 (133) 9.11-15.48 
pump dryer 40 
50 


tv9 
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W vet,plant 


Fires’ ae 


Eg.pp = 
where Ešwh and Exes are the total exergy inputs to the plant at wellhead and reservoir 
conditions, respectively, while they are calculated using Eqs. (6c) and (7). 

Note that the exergy input to the plant is sometimes taken as the exergy of geothermal 
fluid in the reservoir. Those who prefer this approach argue that a realistic and most 
meaningful comparison between geothermal power plants can only be performed 
considering the methods of harvesting the geothermal fluid. However, others argue that 
taking the reservoir as the starting point is not proper for geothermal power plants since 
conventional power plants are evaluated on the basis of the exergy of the fuel burned at the 
plant site [118,119]. 

Exergetic efficiencies of some geothermal power plants worldwide in order of increasing 
efficiency are summarized by the author, as illustrated in Table 10 [119,120-127]. It is 
obvious from this table that exergy efficiencies vary from 16.3% to 53.9% depending on 
the dead state temperature and the technology used. 

DiPippo [120] obtained the following correlation for the exergetic efficiency as a 
function of the dead-state temperature using the Brady geothermal power plant: 


Brady = 23.177 — 0.14297 9. (138) 


DiPippo [120] concluded as follows: “Geothermal binary plants are relatively poor 
converters of heat into work. Thermal efficiencies typically lie in the range of 8-12%. As a 
consequence, a 1-2% point improvement in power output translates into a gain of about 
10-20% in efficiency. The results show that binary plants can operate with very high 
exergetic efficiencies even when the motive fluids are low-temperature and low-exergy. 
Exergetic efficiencies of 40% or greater have been achieved in certain plants with geofluids 
having specific exergies of 200 kJ/kg or lower. The main design feature leading to a high- 
exergy efficiency lies in the design of the heat exchangers to minimize the loss of exergy 
during heat transfer processes. Another important feature that can result in a high exergy 
efficiency is the availability of low-temperature cooling water that allows a once-through 
system for waste heat rejection.” 

In the following, the Grassmann diagram of a 12.4MW existing binary geothermal 
power plant is presented as an illustrative example of geothermal power plants. Exergy 
analysis of this plant was performed by Kanoglu [119] using actual plant data to assess the 
plant performance and the Grassmann shown in Fig. 25 [119] was drawn. As seen in this 
figure, the causes of exergy destruction in the plant included the exergy of the working fluid 
lost in the condenser, the exergy of the brine reinjected, the turbine-pump losses, and the 
preheater—vaporizer losses. The exergy destruction at these sites accounted for 22.6%, 
14.8%, 13.9%, and 13.0% of the total exergy input to the plant, respectively. 

Dagdas et al. [127] proposed a new model in order to improve efficiency and augment 
power output of the single flash type geothermal power plant operating in Denizli- 
Kızıldere, Turkey. In this model, electric energy requirement of the region would be met 
and performance of the geothermal power plant would be up rated by building a gas 
turbine power plant near the existing plant. If the proposed combined cycle power plant 
(gas turbine integrated geothermal power plant) would be operated, power generated in the 
geothermal plant could increase 43.5%. Energy and exergy efficiencies of the plant would 
be increased from 4.56% to 6.5%, and from 19.78% to 28.4%, respectively. 


Table 10 
Geothermal power plant exergetic efficiencies (in order of increasing efficiency) [119,120—-127] 
Technology Plant name Location/country Dead state temperature Exergy efficiency” (%) Sources 
CC) 
Binary Brady bottoming Churchill County, 16.8 16.3 DiPippo [120] 
(p.m. data) NV/USA 
Binary Brady bottoming Churchill County, 30.1 17.9 DiPippo [120] 
(a.m. data) NV/USA 
Binary: recuperated Rotokawa Rotokawa/New Zealand 18.7 DiPippo [120] 
Single flash Kizildere Denizli/Turkey 16 19.78" Dagdas et al. [121,127] 
19.97° 
Single flash Kizildere Denizli/Turkey 16 20.5 Ozturk et al. [122] 
Single flash Kizildere Denizli/Turkey 15 20.8 Cerci [123] 
Single flash Kizildere Denizli/Turkey 18 24.3? Yildirim and Gokcen 
27.2° [124] 
Binary Nigorikawa (Mori) pilot Hakodate on 13 21.6 DiPippo [120] 
Hokkaido/Japan 
Binary Kalina KCS 34, Husavik Husavik/ Iceland 5 23.1 DiPippo [120] 
Single flash Kakkonda Japan 26.5 Setiawan [125], Yildirim 
and Gokcen [124] 
Binary: simple Rotokawa Rotokawa/New Zealand 27.8 DiPippo [120] 
Single flash Otake Japan 29.0 Setiawan [125], Yildirim 
and Gokcen [124] 
Dual-level binary Northern NV/USA 12.8 29.14 Kanoglu [119] 
34.2° 
Single flash Salak (Unit 1&2) Indonesia 30.4 Setiawan [125], Yildirim 


and Gokcen [124] 


99 
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Binary: simple’ 


Single flash 
Single and double flash 
Single and double flash 
Single flash 


Hybrid flash-binary 
Binary: dual-level 


Double flash 
Dry steam 


Dry steam 


Binary: flash evaporator 


Kavala 


Ahuachapan 
(Unit 1 & 2) 
Cerro-Prieto 


Ahuachapan 
(Unit 1 & 2) 
Blundell 


Rotokawa 
Heber SIGC 


Beowawe 
Kamojang 
(Unit 2 & 3) 
Darajat 


Otake pilot 


Eratino/Greece 20 


El Salvador 
Mexico 
El Salvador 


Near Milford in Beaver 
County/USA 

Rotokawa/New Zealand 
Imperial Valley, CA/ 15 
USA 

NV/USA 

Indonesia 


Indonesia 


Island of Kyushu/Japan 18 


Koroneos et al. [126] 


Setiawan [125], Yildirim 
and Gokcen [124] 
Setiawan [125], Yildirim 
and Gokcen [124] 
Setiawan [125], Yildirim 
and Gokcen [124] 
DiPippo [120] 


DiPippo [120] 
DiPippo [120] 


DiPippo [120] 

Setiawan [125], Yildirim 
and Gokcen [124] 
Setiawan [125], Yildirim 
and Gokcen [124] 
DiPippo [120] 


“Exergy efficiencies except those denoted in the above table are based on reservoir exergy. 


At reservoir conditions. 
“At wellhead conditions. 


‘Based on the exergy of the geothermal fluid at the vaporizer inlet. 


“Based on the exergy drop of the brine across the vaporizer-preheater system (i.e., exergy input to the Rankine cycle). 
‘For the exploitation of the field, a Rankine cycle is proposed by Koroneos et al. [126]. 
For the first deposit based on the two case studies. 


For the second deposit. 
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Brine 
reinjection 


Turbine-Pump 
losses 


Vaporizer- 
Preheater 
losses 


Turbine 844 kW, 13.9% 


Pump 15 kW, 0.25% Condenser 


losses 


Brine 
exergy input 
6073 kW 


1374 kW 
22.6% 


390 kW Parasitic power 
64% =—————— 


1769 kW 
29.1% Net power 


Fig. 25. Grassmann (exergy flow) diagram of a 12.4MW existing binary geothermal power plant drawn by 
Kanoglu [119]. 


4.4. Biomass 


Chemical exergy of a substance is the maximal possible useful work that may be 
produced by process of physical and chemical equilibration of the substance with the 
ambient. To evaluate the chemical exergy of a mixture, the knowledge of its enthalpy of 
combustion, elemental composition and absolute entropy is necessary. When some of these 
data are absent, the methods of estimation of energy are used [128]. The chemical exergy of 
gaseous and some liquids can be obtained using standard chemical exergy tables given in 
Ref. [7] if the chemical composition of the fuel is known. Many technical liquid and solid 
fuels are multi-component solutions (or mixtures) of very complicated, usually unknown 
components. An exact solution of the chemical exergy of such fuels is impossible. In this 
regard, Rant’s relation including a constant ratio of chemical exergy and calorific value (or 
heating value) for liquid and (separately) for solid fuels (Excy/HV) may be used. 
However, the calculation for different organic substances indicated that the ratio 
(Excy/HV) depends significantly on the chemical composition. 

If the elemental composition is known, the chemical exergy of various fuels can be 
evaluated accurately using the Szargut’s correlations listed elsewhere [7]. In addition, a 
correlation formula to estimate chemical exergies of oil fractions and fuel mixtures from 
enthalpy of combustion and atomic composition developed by Govin et al. [128] is used. 

Based approximation method for fuel chemical exergy (Yen), the following relations are 
widely used. 


Excy 


Buy = pa (139) 
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Excy 
Bunny = HHV’ 
where f is the proportionality constant (or quality factor or exergy coefficient), while LHV 
and HHV denote the lower heating (net caloric) and higher heating (gross caloric) values, 
respectively [129]. 
Szargut [22] proposed that rnv and Byyy may be taken to be approximately 1.15 and 
1.05 for wood, respectively. The following relations for wood with the limit of 
(Zo, /Zc)<2.67 are proposed by Szargut et al. [7] and Szargut [22], respectively. 


Buav = {1.0412 + 0.2160(Zu,/Zc) — 0.2499(Zo,/Zc)[1 + 0.7884(Zn, /Zc)] 
+0.0450(Zn, /Zc)}/{1 — 0.3035(Zo,/Zc)} (141) 


(140) 


and 


Piny = {1.0412 + 0.2160(Z1, /Zc) + 0.2499(Zo,/Zc)[1 + 0.7884Zn, /Zc)] 
+0.0450(Zn,/Zc)}/{1 + 0.3035(Zo,/Zc)}. (142) 


where Z; is mass fraction of element i. 

Eq. (141) has been successfully used by some investigators [130-132]. In this regard, 
Nilsson [130] used the following relation for determining the specific chemical exergy of 
straw 


EXcH straw = (LHV F LwZw) Buay F ExcuwZw, (143) 


where L and Z denote the enthalpy of phase change (kJ/kg) and the mass fractions (—), 
respectively. Using the following data for moist straw; Zc = 38.9%, Zp, = 4.9%, 
Zo, = 33.8%, Zn, = 0.4%, Zy = 18.0%, LHV = 13,800 kJ/kg, Ly = 2440kJ/kg and 
Excu.w = 50kJ/kg, the specific chemical exergy was calculated to be 16,100kJkg. In 
other words, the ratio of chemical exergy to the lower heating value was 1.16. In the 
calculation, the reference temperature was 25°C, the relative humidity of atmospheric air 
was 70% and the chemical exergy of ash was neglected. 

Panopoulos et al. [132] analyzed exergetically the combination of allothermal biomass 
gasification and solid oxide fuel cell for small scale CHP, and presented exergy efficiency 
relations. They also determined that for SOFC operation at current density 2500 A/m”, the 
system used 90kg/h biomass and operated with electrical exergetic efficiency 32% 
producing 140kW,, while the combined electrical and thermal exergetic efficiency was 
35%. 

Franco and Giannini [133] analyzed the perspective of improving the efficiency of 
biomass for power generation and proposed an exergy loss based economic analysis of 
biomass utilization. A thermoeconomic method of analysis based on the consideration of 
the cost of the exergy losses showed how the conversion efficiency of renewable sources as 
biomass, had a lower limit that, in a lot of cases, was considerably higher than the actual 
efficiency level. 

According to Shieh and Fan [134], specific chemical exergy (kJ/kg) of structurally 
complicated materials, e.g. biomass, can be estimated from their elemental compositions as 
given below [135] 


Won = 4.1868{8177.79[C] + 5.25[N] + 27, 892.63[H] — 3173.66[0] 
+ 0.15{O](7837.677[C] + 33, 888.889[H] — 4236. 1[O]}. (144) 
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Table 11 
The calculated exergy values of various of woody biomass types in Turkey [136] 


Type of wood Ultimate analysis (%) Exergy (MJ/kg) HHV (MJ/kg) Buy 


C H O N 
Hus wood 48.5 5.6 45.7 0.2 17.90 17.22 1.040 
Willow wood 50.8 6.0 42.5 0.6 19.63 19.01 1.033 
Sugar cane 47.3 6.2 46.2 0.3 18.15 17.52 1.036 
Olive residue 49.5 6.3 43.7 0.5 19.39 18.77 1.033 
Straw 47.0 6.1 46.4 0.5 17.90 17.25 1.038 
Corn cob 46.2 6.3 47.2 0.3 17.75 17.12 1.037 
Olive shell 52.4 6.0 40.4 1.5 20.45 19.76 1.031 
Wheat straw 48.7 6.2 44.8 0.3 18.83 18.21 1.034 
Rice shall 48.7 6.0 44.7 0.6 18.60 17.96 1.036 
Rice straw 47.0 6.4 45.5 1.1 18.38 17.77 1.034 
Hybride poplar 51.6 6.3 41.5 0.6 20.42 19.84 1.029 


Cotton stalk 46.7 6.2 46.5 0.6 17.90 17.96 1.037 


Bilgen et al. [136] studied on exergy analyses of various types of wood and coal samples 
with Turkish origins by using exergy analysis method developed and published by 
Srivastava [137] and Bejan et al. [138]. They calculated the exergy values of various of 
woody biomass types in Turkey, as given in Table 11 [136]. As can be seen in this table, the 
values for Byyy vary from 1.029 to 1.040. On the other hand, these values were found 
greater than 1.0 for all the woody biomass types. Bilgen et al. [136] concluded that this 
result proposed that maximum work is always greater than the heating value, or energy 
delivered by a process in which these coals are used as a functional material. 

Utlu et al. [139] evaluated the Turkish cotton stalk production in 2003 using energy and 
exergy analysis method for the first time. The values for the net energy and exergy gained 
were obtained to be about 49,146 and 38,426 MJ/ha, respectively. The overall mean energy 
and exergy efficiencies of the cotton production in the year studied were found to be 
33.06% and 33.12%, respectively. 

Prins et al. [140] performed laboratory experiments in a small-scale (5-10 g sample) 
fixed-bed torrefaction reactor placed inside an oven and gave exergetic efficiency 
relations for various processes. They reported as follows: Overall exergetic efficiency of 
air-blown gasification of torrefied wood was lower than that of wood since the volatiles 
produced in the torrefaction step were not utilized. For the entrained flow gasifier, the 
volatiles could be introduced into the hot product gas stream as a ‘chemical quench’. The 
overall efficiency of such a process scheme was comparable to direct gasification of wood, 
but more exergy was conserved in as chemical exergy in the product gas (72.6% versus 
68.6%). 


4.5. Other renewable energy systems 


The following section summarizes exergy relations for ocean surface waves, pre- 
cipitation, ocean thermal gradient and tides, which have been taken from a study 
conducted by Hermann [14]. These relations may be used in exergetic evaluating these 
types of RERs. 
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4.5.1. Ocean surface waves 

The average specific exergy of surface waves ranges from 10 to 100kW/m, with much 
higher values during storms [14,141]. The waves concentrate kinetic energy transfer from 
the wind over many kilometers and can deposit this energy at the coastline. The primary 
factors determining the exergy of a fluid gravity wave are the period and [14,142]. For an 
approximately sinusoidal deep ocean wave, the length-specific exergy flux is 


1 2 
Wwave = 7 PV waveIGvaves (145) 
where Vwave is the wave speed, g is the gravitational constant and wave is the amplitude. 


4.5.2. Precipitation 

The specific exergy of precipitation can be broken down into physical and chemical 
components due to its gravitational and diffusive chemical potential energy relative to 
seawater. 

For a specific replenishing body of water such as a dammed reservoir, the gravitational 
potential exergy is 


Y pre,pot = 192, (146) 


where m is the mass flow rate and z is the vertical distance from the water level of the 
reservoir to the reference height or average sea level [14]. 

Precipitation also contains a diffusive exergy component due to the difference in 
chemical potential between the Earth’s freshwater and seawater. Since the difference 
between seawater and freshwater is primarily dissolved salts, this value can be computed 
most directly by considering a reversible osmosis process. 

The work done against the osmotic pressure is 


W pre,diff = —RT Ind — Ysolute)» (147) 


where y is the mole fraction of solute in the water [14,143]. 


4.5.3. Ocean thermal gradient 

The solar warming of seawater in the tropics and the radiative loss of thermal energy to 
space near the poles creates a density-driven current which circulates seawater from the 
surface to the deep ocean. Due to the small temperature differential, the exergy density of 
the ocean thermal gradient is relatively small. For a temperature difference of 20 K, the 
exergy content is approximately 7% of the thermal energy [14]. 


YTEG = mC p,ave(T warm E Teola) (1 _ Toa ) . (148) 
warm 

4.5.4. Tides 

The specific tidal exergy of a designated area of ocean or enclosed reservoir is a function 
of the record of minima and maxima and the time periods between these levels [14]. As the 
tide rises and falls, an average exergy flow considering the reservoir area integrated over 
the height difference between local maxima and minima and the period between these 
points can be written as follows: 


Leg fe" Ares(z) dz 
Via a ; (149) 


n 
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where f¢ is the period between local maxima and minima of the tidal record, Ares the 
area of the reservoir at height z, and n the index number for each maximum or 
minimum [14]. 


4.6. Country based renewable energy sources 


The method of exergy analysis has been applied to a wide variety of thermal and thermo 
chemical systems. A particular thermo dynamical system is the society, for example, of a 
country or a region [144]. In order to attain energy saving, the application of energy and 
exergy analysis modeling techniques for energy-utilization assessments has recently 
significantly increased. Utlu and Hepbasli [145] conducted a comprehensive review on the 
exergetic assessment of countries. In this regard, the approaches used to analyze energy 
utilization of countries or societies may be grouped into three types, namely Reistad’s 
approach, Wall’s approach and Sciubba’s approach. In Reistad’s approach, the end use is 
divided into three sectors, these names are industry, transportation and residential-com- 
mercial. In addition, energy sector with electric generation and distribution, oil refining, is 
treated separately. Flows of energy carriers for energy uses are not included in this 
approach. In Wall’s approach, all types of energy and material flows are considered. 
However, to energy carriers for energy use are considered, these flows encompass wood for 
construction materials and for the pulp and paper industry, harvested wood and fodder, 
and the products from these materials. The so-called Sciubba’s approach is relatively a new 
one. Extended-exergy accounting (EEA) method introduced by Sciubba was applied to the 
Italian society in 1996 by Milia and Sciubba [146] and recently to Norway with the figures 
of 2000 by Ertesvag [147]. 

Table 12 [145,149,150] illustrates Eqs. (150)-(158) used to estimate exergy (second law) 
utilization efficiencies for the RERs. The exergy efficiency frequently gives a finer 
understanding of performance than (n). The parameter (e) weights energy flows by 
accounting for each in terms of availability. It stresses that both losses and internal 
irreversibilities need to be dealt with to improve performance. 

To date, various studies, as reviewed in detail elsewhere [145], were performed to 
analyze energy and exergy utilization efficiencies of countries in the four subsectors, 
namely utility, industrial, residential- commercial and transportation sectors. However, 
studies conducted on the evaluation of energy utilization efficiency of RERs using energy 
and exergy analyses are very limited in numbers in the open literature. The works done by 
Hepbasli and Utlu [15,148] on Turkey’s RERs for the years of 2001 and 2010 appear to be 
the first ones. 

Table 13 [15,148,151—-154] illustrates the main outcomes of the exergy utilization analysis 
for both Turkey’s general and RERs. It is clear from this table that the exergy utilization 
efficiencies for Turkey’s general ranged from 13.10 to 27.94%. For Turkey’s RERs, exergy 
utilization efficiencies varied from 24.14 and 44.30% in 2001 and 2010, respectively [15]. 
These results have clearly indicated the necessity of the planned studies towards increasing 
exergy utilization efficiencies in the sector studied. Hepbasli and Utlu [15] concluded that the 
analyses reported here would provide the investigators with knowledge about how much 
effective and efficient a country uses its natural resources, while this knowledge is also needed 
for identifying energy efficiency and/or energy conservation opportunities as well as for 
dictating the energy strategies of a country or a society, especially for a sustainable future. 


Table 12 
Relations used in the exergy analysis [145,149,150] 
Equations Description 
1 e = (exergy in products/total exergy input) x 100 (150) Definition of exergy utilization efficiencies 
2 e = (n/B)1 — [To /(T — To) In(T/To)} (151) Exergy utilization efficiencies of space heating 
3 e = (Wou/Win)100 (152) Exergy utilization efficiencies of RERs for utility sector 
4 e = nil — (To/ T) (153) Exergy utilization efficiencies of cooking 
5 e=m(To/T) - 1] (154) Exergy utilization efficiencies of refrigeration 
6 foe = [(a1 £a1) + (a28a2) + (a3£a3) +--+ + (a7£a7)]/100 (155) Overall sector exergy utilization efficiencies for electric use 
7 for = [Pan X ent) + Pon X ewhe) + Fe X &,1)]/100 (156) Overall sector exergy utilization efficiencies for fuel use 
8 Eor = [("sh X €2sh,or) + (wh X Ewh,orf) + (Te X &2e,r0)]/('sh + wn +r) (157) Overall sector exergy utilization efficiencies for RE use 
9 £o = (ĉc€rc + £off exre + Erlexre)/(Cre + f exre + exre) (158) Overall sector exergy utilization efficiencies for residential-commercial sector 
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Table 13 
Comparison of total exergy inputs/outputs as well as exergy utilization efficiencies of Turkey’s renewables and total [15,148,151—154] 
Investigators Year Population Total exergy input Total exergy output RERs* input? RERs* output Efficiencies 
analyzed/ 
published 
(PJ) (GJ/cap) (PJ) (GJ/cap) (PJ) (GJ/cap) (PJ) (GJ/cap) Total RERs*° 
e (%) €(%) 
Unal [151] 1991/1994 57,024,515 2279 39.97 539.5 9.46 450.54 7.90 23.7 
Rosen and Dincer [152] 1993/1997 58,808,625 680.6 11.57 445.2 7.57 455.1 7.74 27.1 
Ileri and Gurer [153] 1995/1998 59,706,545 2697.3 45.17 352.3 5.9 4714.4 795 13.1” 
Utlu and Hepbasli [154] 1999/ 66,022,636 3380.34 51.2 499.61 7.55 459.25 6.96 24.04 
2000/ 67,803,927 3469.62 52.55 525.77 7.96 449.72 6.63 24.78 
Hepbasli and Utlu [15] 2001/2004 68,820,985 3139.07 45.61 783.75 11.39 395.68 5.75 95.51 1.38 25.22 24.14 
Hepbasli and Utlu [148] 2010/ 76,840,418 7400.34 96.31 2067.65 26.91 580.21 7.55 351.7 4.58 27.94 44.30 


“These values are found by calculation from references. 
PExcluding utility sector use. These values are obtained to be 43.61 (21.83)% with the utility sector use of 236.6 PJ. 


“Renewable energy resources. 
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5. Conclusions 


Exergetic aspects of RERs are comprehensively presented in this study. These RERs 
studied are solar, wind and geothermal energy systems as well as biomass and country 
based RERs. Studies conducted on these RERs are also compared with the previously ones 
in tabulated forms. 

Some concluding remarks which can be extracted from this study are as follows: 


(a) Exergy is a way to a sustainable development. In this regard, exergy analysis is a very 
useful tool, which can be successfully used in the performance evaluation of RERs as 
well as all energy-related systems. 

(b) Based on the review results of the RERs, the following exergy efficiency values were 
found. 

@ The exergy efficiency of a SPC was found to range from 0.4% to 1.25% [33]. 

e The exergy efficiencies of a solar collector, a PV and a hybrid solar collector 
were found to be 4.4%, 11.2% and 13.3%, respectively. Solar collector had 
the lowest value. In terms of energy quality, the hybrid solar collector was very 
effective [67]. 

e Based on an exergetic evaluation of wind energy systems, average differences 
between energy and exergy efficiencies were obtained to be approximately 40% at 
low wind speeds and up to approximately 55% at high wind speeds [68]. 

e@ Geothermal resources should be classified by their exergy, a measure of their ability 
to do work. The classification with reference to SExI was more meaningful as there 
was no general agreement on the arbitrary temperature ranges used in the 
classification of geothermal resources by temperature [70,71]. 

e Based on a comparison of exergy efficiency values for GDHSs installed in Turkey, 
exergy efficiencies ranged from 42.94% to 66.13% depending on the various dead 
state temperatures. As expected, the lower the reference state temperature, the 
significantly larger the exergy losses in the system pipeline and heat exchangers 
[8 1-83,86,87,90,91]. 

e Exergy efficiency values of GSHP units varied between about 67% and 92% on 
the product/fuel basis depending on the various dead state temperatures 
[93-97,103,106,107]. 

e It was found from an exergetic evaluation of a drying process utilizing a GSHP 
system that the exergy efficiency values of the drying process ranged from 81.35% to 
87.48% based on the inflow, outflow and loss of exergy, and 9.11% to 15.48% based 
on the product/fuel basis between the same drying air temperatures with a drying air 
mass flow rate of 0.12 kg/s. It was also concluded that it was more meaningful, 
objective and useful to assess the performance of the drying process relative to the 
performance of similar drying processes on the product (or benefit)/fuel basis [115]. 

e Exergetic efficiencies of some geothermal power plants worldwide were compared 
with each other and found to vary from 16.3% to 53.9% depending on the dead 
state temperature and the technology used [119,120,122—127]. 

e Exergy utilization efficiencies of Turkey’s RERs were found to range from 24.14% 
and 44.30% in 2001 and 2010, respectively [15]. These results have clearly indicated 
the necessity of the planned studies towards increasing exergy utilization efficiencies 
in the sector studied. 
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(c) As a conclusion, the author expects that the analyses reported here will provide the 
investigators with knowledge about how effective and efficient a country uses its RESs. 
This very useful knowledge is also needed for identifying energy efficiency and/or 
energy conservation opportunities, as well as for dictating the right energy and exergy 
management strategies of a country. 
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